AdOO  31U 


COMPUTING  INTERNAL  COCKPIT  REFLECTIONS  OF  EXTERNAL 
POINT  LIGHT  SOURCES  FOR  THE  MODEL  YAH-64  ADVANCED 
ATTACK  HELICOPTER  (LOW  GLARE  CANOPY  DESIGN) 


Christopher  C.  Smyth 


July  1977 

AMCMS  Code  624209 .C5 204 12 


Approved  for  public  release, 
distribution  unlimited. 


U ! 


U,  S.  ARMY  HUMAN  ENGINEERING  LABORATORY 
Aberdeen  Proving  Ground,  Maryland 


Destroy  this  report  when  no  longer  needed. 

Do  not  return  it  to  the  originator. 

The  findings  in  this  report  are  not  to  be  construed  as  an  official  Department 
of  the  Army  position  unless  so  designated  by  other  authorized  documents. 

Use  of  trade  names  in  this  report  does  not  constitute  an  official  endorsement 
or  approval  of  the  use  of  such  commercial  products. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Wian  Data  Bnlarad) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER 

Technical  Memorandum  24-77 

2.  GOVT  ACCESSION  NO 

3.yrt^lPIENT‘S  CATALOG  NUMBER 



- ■».-  TITLE  '■  

COMPUTING  INTERNAL  COCKPIT  REFLECTIONS  OF 
;^XTERNAL  POINT  LIGFTT  SOURCES  FOR  THE  MODEL 
XAH-64  advanced  attack  HELICOPTER  L^W 

.GCARE  Canopy  design),  ' 

i5'-TTPCTfrRP^ttT'^”PLHII^iyeRED 

4.  ncnwBRMinJTrrroi  report  nLmber 

Christopher  C^myth  j 

B.  CONTRACT  OR  GRANT  NUMBERC*; 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

US  Army  Human  Engineering  Laboratory 
Aberdeen  Proving  Ground,  Maryland  21005  ^ 

10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  & WORK  UNIT  NUMBERS  ' 

AMCMS  Code  624209 .C52041 2 

* 

M.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

(4 

61 

U.  MONITORING  AGENCY  NAME  8 AOORESS^I/  dltUtmt  from  ControWng  OUleo) 

IS.  SECURITY  CLASS,  (ol  Ihit  report) 

Unclassified 

IS*.  DECLASSIFICATION/ downgrading 
SCHEDULE 

t«.  distribution  statement  Cat  t/ila  Raport) 


Approved  for  public  release:  distribution  unlimitnd. 

'f'm 

17.  OISTRIBUTW#- 


t 


im  iriinniffiSnSxirTT^  20,  tt  dJtUrlnt  from  R0pi,r1) 


18.  supplementary  notes 


19.  KEY  WORDS  (Continu*  on  r«v«r«>«  It  n«c««««ry  «nd  ld*n(l/y  by  block  numbor) 


Internal  Cockpit  Reflections 
Primary  Reflections 
External  Point  Light  fiources 
Transparent  Canopy  Surfaces 
Canopy  Surfaces 


Canopy  Layouts 
Helicopters 

Model  YAH-64  Advanced  Attack  Helicopter 
Human  Factors  Engineering 


ABSTRACT  (Voatfauo  on  oM*  tt  nmemnmarr  md  Idotutty  by  block  numbor) 


The  US  Army  Human  Engineering  Laboratoy  (HEL)  has  developed  a computer  program  for 
computing  the  inteina!  cockpit  reflections  on  the  transparent  canopy  surfaces  of  external  point 
light  sources.  Computations  have  been  completed  for  the  Model  YAH-64  Advanced  Attack 
Helicopter  (low  glare  canopy  design).  The  results  show  that  primary  reflections  as  seen  from  the 
pilot's  position  are  possib'e  on  (1 ) the  upper  rear  corners  of  the  forward  side  canopy  surfaces,  (2) 
the  upper  edges  of  the  rear  sides,  and  (3)  the  sides  of  the  top  surface.  Computations  have  also 
been  completed  for  the  copilot’s  position  and  show  possible  reflections  on  the  front  and  side 
surfaces.  A computer  graphics  output  is  used  to  show  reflection  points  on  canopy  layouts  and 
_ perspectives  of  the  cockpit. 

DO  /jTn  U73  EDITION  OF  I 65  IS  OBSOLETE 


. " SECURITY  CLASSIFICATION  OF  THIS  PAGE  (WhOTi  D«l«  Bntofd) 


AMCMS  Code  624209.C52041 2 


Technical  Memorandum  24-77 


COMPUTING  INTERNAL  ^vOCK^IT  REFLECTIONS  OF  EXTERNAL 
POINT  LIGHT  SOURCES  FOR  THE  MODEL  YAH-64  ADVANCED 
ATTACK  HELICOPTER  (LOW  GLARE  CaNOPY  DESIGN) 


Christopher  C.  Smyth 


July  1977 


US  ARMY  HUMAN  ENGINEERING  LABORATORY 
Aberdeen  Proving  Ground,  Maryland  21005 


I 

I 

I 


Approved  tor  public  release, 
distribution  unlimileit 


1 


CONTENTS 


INTRODUCTION  . . 3 

METHOD  3 

DISCUSSION  5 

FURTHER  RESEARCH 19 

CONCLUSION  19 

REFERENCES 26 

APPENDIXES 

A.  Internal  Reflections  From  Cylindrical  Surfaces 27 

B,  Computer  Program 33 

FIGURES 


1 . Top,  Side  and  Front  Views  of  Canopy  Frame  and  Obstructing  Surfaces  6 

2.  Top,  Side  and  Front  Views  of  Canopy  Frame  7 

3.  Perspective  View  of  the  Cockpit  Interior  From  the  Pilot’s  Nominal 

Viewing  Position  and  Direction 8 

4.  Entry  Ray  Positions  Generating  Primary  Reflections  on  the  Right-Hand 

Side  of  the  Canopy  As  Seen  From  the  Pilot’s  Position  . 9 

5.  Primary  Reflection  Points  on  the  Right-Hand  Side  of  the  Canopy  and 

Their  Associated  Reflectance  Values  for  the  Pilot’s  Position 10 

6.  Perspective  View  of  Entry  Ray  Positions  for  the  Pilot’s  Nominal  Viewing 

I Direction  and  Position  11 

' 7.  Perspective  View  of  Primary  Reflection  Points  for  the  Pilot’s  Nominal 

Viewing  Direction  and  Position 12 

8.  Perspective  View  of  Primary  Reflection  Points  for  the  Pilot  Viewing 

20-Degrees  to  the  Right  Side 13 

9.  Entry  Ray  Positions  on  the  Right  Rear  Side  Canopy  Surface  and  Their 

Corresponding  Reflection  Points  on  the  Right  Side  of  the  Top  Surface 14 

10.  Entry  Ray  Positions  on  the  Right  Forward  Side  Canopy  Surface  and  Their 

Corresponding  Reflection  Points  on  the  Right  Front  of  the  Top  Surface  .....  15 

1 1 . Entry  Ray  Positions  on  the  Left  Rear  Side  Canopy  Surface  and  Their  Corresponding 

Reflection  Points  on  the  Top  Edge  of  the  Right  Rear  Surface  16 

12.  Entry  Ray  Positions  on  the  Lower  Front  Canopy  Surface  and  Their  Corresponding 

Reflection  Points  on  the  Upper  Rear  Corner  of  the  Right  Forward  Side  Surface  . 1 

13.  Entry  Ray  Positions  on  the  Left  Forward  Side  Canopy  Surface  and  Their  Corresponding 
Reflection  Points  on  the  Upper  Front  Edge  of  the  Right  Rear  Side  Surface  ...  18 

14.  Entry  Ray  Positions  Generating  Reflections  on  the  Right  Hand  Side  of  the  Canopy 

as  Seen  From  the  Copilot’s  Nominal  Position 20 


1 


15.  Primary  Reflection  Points  on  the  Right  Hand  Side  of  the  Canopy  and  Their  Associated 


Reflectance  Values  for  the  Copilot’s  Position 21 

16.  Perspective  View  of  Primary  Reflection  Points  for  the  Copilot’s  Viewing 

45-Degrees  to  the  Right  Side 22 

1 7.  Entry  Ray  Positions  on  the  Right  Forward  Side  Canopy  Surface  and  Their  Corresponding 

Reflection  Points  on  the  Right  Side  of  the  Upper  Front  Surface  23 

18.  Entry  Ray  Positions  on  the  Right  Forward  Side  Surface  and  Their  Corresponding 

Reflection  Points  on  the  Upper  Right  Cornet  of  the  Lower  Front  Surface  ....  24 

19.  Entry  Ray  Positions  on  the  Left  Forward  Side  Surface  and  Their  Corresponding 

Reflection  Points  on  the  Upper  Edge  of  the  Right  Forward  Side  25 


2 


COMPUTING  INTERNAL  COCKPIT  REFLECTIONS  OF  EXTERNAL 
POINT  LIGHT  SOURCES  FOR  THE  MODEL  YAH-64  ADVANCED 
ATTACK  HELICOPTER  {lOW  GLARE  CANOPY  DESIGN) 


INTRODUCTION 

The  US  Army  Human  Engineering  Laboratory  (HEL)  has  developed  a compu'er  program  for 
computing  internal  cockpit  reflections  on  the  transparent  canopy  surfaces  of  external  point  light 
sources.  This  work  is  part  of  a three-stage  effort  to  determine  optimum  canopy  designs  for  the 
Model  209  AH-1S  Cobra  Helicopter  and  the  Model  YAH-64  Advanced  Attack  Helicopter  (AAH). 
This  work  was  undertaken  at  the  request  of  the  Project  Manager’s  Office,  USA  Aircraft 
Survivability  Equipment.  The  low  glare  canopy  design  presently  used  on  both  models  consists  of 
flat,  transparent  panels  on  the  front  surfaces  and  simple  cylindrical  panels  on  the  sides  and  top. 
The  design  is  a reasonable  choice  for  reducing  both  solar  glint  to  outside  observers  during 
daytime  operations  and  internal  reflections  of  outside  light  sources  during  nighttime  operations. 

A flat  plate  canopy  (FPC)  design  was  originally  developed  for  the  Cobra  and  AAH  to  reduce 
dayt'me  solar  glint  to  a momentary  flash  at  certain  obseiver-aircraft  sun  angles.  A moving 
aircraft  no  longer  produced  the  continual  solar  glint  which  was  present  on  the  earlier 
compound-shaped  canopy  designs.  The  continual  presence  of  solar  glint  had  increased  the  range 
of  visual  detection  by  ground  observers. 

However,  in  certain  lighting  situations  during  nighttime  operations,  the  internal  surfaces  of 
the  FPC  performed  as  mirrors  reflecting  virtual  images  of  external  Pght  sources  that  were  visible 
to  the  pilot.  HEL  has  shown  by  computer  analysis  that  these  refic  .tioris  are  possible  on  most  of 
the  transparent  surfaces  and  for  a wide  range  of  source  locations  (5).  The-se  virtual  images  of 
ground-level  lights  were  disorienting  to  the  pilot  and  he  could  not  easily  discriminate  between  the 
light  sources  on  the  ground  and  their  reflections  from  the  canopy  surfaces.  This  problem  was  a 
potential  safety  hazard  during  flight. 

The  present  low  glare  canopy  design  was  developed  to  reduce  these  two  conflicting 
problems  to  manageable  levels.  The  design  incorporates  front  planar  transparent  surfaces  and 
simple  cylindrical  surfaces  for  the  sides  and  top.  HEL  recommended  a similar  design  with, 
however,  novel  features  (6).  The  present  work  effort  is  directed  toward  a closer  study  of  the  two 
problems  of  glint  and  reflections,  and  developing  an  optimum  design  for  the  canopy’s  transparent 
surfaces. 


METHOD 

A ray-tracing  program  was  written  to  trace  in  three  dimensions  the  straight-line  rays  from 
the  nominal  position  of  the  pilot’s  eye  backwards  to  visible  points  on  the  internal  surfaces  of  the 
cockpit.  Each  ray  is  traced  between  transparent  surface  points  until  a nontransparent  surface  is 
reached.  These  surfaces  are  assumed  to  be  diffusive  without  specular  reflectances  and  the  ray  is 
considered  absorbed.  At  each  reflection  point  on  a transparent  surface,  the  reflectar  :e  and 
transmittance  are  computed  along  with  the  directional  cosines  of  the  corresponding  transmitted 
and  reflected  rays.  In  this  way,  a reflected  ray  reaching  the  pilot’s  eye  is  traced  backwards  to  all 
possible  externa!  sources  that  can  generate  that  ray. 
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The  transparent  surfaces  of  the  low  glare  canopy  design  are  specified  as  a set  of  planar  and 
cylindrical  surfaces  and  their  corresponding  edge  vertices.  Each  planar  surface  is  specified  by  the 
coordinates  of  its  edge  vertices  and  the  consecutive  order  in  which  adjacent  vertices  are  listed.  A 
cylindrical  surface  is  specified  by  cylindrical  parameters  and  the  consecutive  seauence  of  the  edge 
vertices  and  their  coordinates.  The  cylindrical  parameters  are  (1)  origin  point  on  the  cylindrical 
axis,  (2)  directional  cosines  of  the  axis,  and  (3)  the  radius  of  the  cylinder.  The  edges  of  the 
cylindrical  surface  are  assumed  for  simplicity  to  be  curvilinear  lines  which  become  straightened 
when  the  cylinder  is  transformed  into  a flat  plane. 

Given  dirtctional  cosines  and  an  origin  point  of  a straight-line  ray,  the  program  computes  in 
turn,  the  intersection  point  of  the  ray  with  each  surface.  The  program  tests  the  intersection  point 
against  the  surface  edges.  The  reflection  point  for  the  ray  is  that  intersection  point  which  is 
contained  within  the  edges  of  the  corresponding  surface  segment.  The  angle  of  incidence  between 
the  surface  normal  and  the  ray  at  this  point  is  computed  along  with  the  corresponding  values  of 
reflectance  and  transmittance  and  the  directional  cosines  of  the  transmitted  and  reflected  rays. 
Tracing  backwards,  the  reflected  ray  becomes  the  incident  ray  for  the  next  set  of  computations. 
(See  Appendix  A for  derivation  of  equations  used  in  ray  tracing  on  cylindrical  surfaces  and 
Appendix  A of  reference  5 for  ray  tracing  on  planar  surfaces  and  computation  of  transmittance 
and  reflectance  values.) 

The  program  includes  interna!  and  external  obstructing  surfaces  and  the  internal  blast 
barrier  of  the  YAH-64  between  the  pilot  and  copilot,  as  well  as  the  transparent  surfaces  of  the 
canopy  in  the  computation.  The  obstructing  surfaces  are  those  that  either  obstruct  the  pilot’s 
vision  or  block  incident  rays  from  external  sources.  The  internal  surfaces  are  (1)  the  pilot’s  seat, 
display  panel  and  side  armor,  (2)  the  copilot’s  seat,  gunner-sight  and  side  armor,  and  (3)  the  sides 
and  floor  of  the  cockpit.  The  external  surfaces  are  (1)  aircraft  nose  section,  (2)  gun  pods  and 
wheel  wells,  (3)  wing  stubs,  (4)  rocket  pods,  (5)  engine  intakes,  and  (6)  rotary  housing  These 
surfaces  are  specified  as  planar  segments  in  the  same  manner  as  are  the  canopy  surfaces.  The 
intersection  computations  are  performed  first  for  all  obstructing  surfaces  and  computation  of  a 
reflection  point  f^or  a ray  on  an  obstructing  surface  renders  the  computation  complete  since  the 
backwards  traced  ray  is  considered  absorbed. 

The  transparent  blast  barrier  which  separates  the  copilot  and  pilot,  is  treated  first  as  a 
reflecting  surface  and  then  as  a transmitting  surface  for  reflection  points  on  surfaces  beyond  it. 

This  computation  process  is  repeated  for  pilot-viewing  directions  indexed  at  equal 
increments  over  a quarter  sector.  The  sector  is  bounded  by  vision  directly  to  the  front,  to  the 
side,  top  and  bottom.  In  this  way,  a table  is  constructed  which  lists  at  discrete  intervals  all 
possible  internal  reflection  points  and  the  corresponding  external  light  directions.  This  approach 
generates  a large  amount  of  data  and  a computer-graphics  routine  is  included  for  output.  The 
primary  reflection  points  and  the  corresponding  incident  ray  entry  points  are  shown  on  side,  top 
and  front  views  of  the  canopy  and  on  perspective  drawings  of  the  cockpit  as  seen  from  the  pilot’s 
position.  Similar  comments  apply  to  computations  for  the  cop  lot’s  position.  (See  Appendix  B of 
reference  5 for  a discussion  of  perspective  drawings.) 
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DISCUSSION 


The  results  of  this  application  are  shos^n  in  Figures  1 through  19.  These  figures  are  hard 
copies  of  the  computer  graphics  output.  Figure  1 shows  side,  top  and  front  viev*  • of  the  canopy 
frame,  blast  barrier  and  obstructing  surfaces.  The  pilot’s  nominal  eye  position  i shown  in  each 
view.  The  blast  barrier  and  obstructing  surfaces  are  sketched  in  with  broken  lin  is  The  aircraft 
fuselage  and  tail  assembly  are  not  included  in  this  sketch. 

Figure  2 shows  side,  top  and  front  views  of  the  canopy  frame  and  blast  barrier,  sketched 
with  broken  lines,  separated  from  the  obstructing  surfaces.  Figure  3 is  a perspective  drawing  of 
the  cockpit  as  seen  from  the  pilot’s  position.  The  pilot’s  nominal  viewing  direction  is  shown  by 
the  small  cross  near  the  top  center  of  the  upper  front  canopy  surface.  The  drawing  covers  a 
60-degree  field  of  view  and  shows  that  the  lower  portions  of  the  front  and  forward  side  canopy 
surfaces  are  blocked  from  view  by  the  pilot’s  instrument  panel. 

The  frame  edges  for  the  canopy  sides  are  drawn  as  straight  lines  connecting  adjacent  corner 
vertices.  This  is  done  for  convenience  in  the  computer  graphics  routines.  The  computations 
assume  that  the  ftame  edges  for  the  cylindrical  surfaces  are  curvilinear  lines  (see  Method,  page  3). 

Figures  4 through  13  show  ' for  the  entry  positions  of  external  rays  generating 
primary  reflections  on  the  right-hand  side  of  the  canopy  for  the  pilot’s  position.  .h'wn  arc 
the  corresponding  primary  reflections  spaced  at  two  degrees  by  two  degrees  increments.  The 
number  shown  at  each  reflection  point  is  equal  to  the  negative  value  of  toe  logarithm  (base  10)  of 
the  I'ght  reflectance.  The  numbers  are  truncated  to  their  integer  values  by  dropping  the 
fractional  parts.  The  numerical  “zero”  corresponds  to  those  reflectaixes  which  are  greater  than 
0.1  in  value.  The  numerical  “one”  corresponds  to  those  values  equal  to  or  le.s  than  0.1  but 
greater  than  0.01 . 


Figure  4 shows  that  entry  points  are  possible  over  much  of  the  lower  front  panel  and  the 
side  surfaces.  Figure  5 shows  that  primary  reflection  points  can  occur  on  (1)  the  upper  rear 
corner  of  the  front  side  panels,  (2)  the  upper  edge  c*'  the  rear  side  panels,  and  (3)  the  side  edges 
of  the  top  panel.  The  front  side  panel  reflections  have  reflectance  values  in  the  0.1  to  1.0  range, 
while  those  on  the  rear  sides  and  top  are  in  the  0.01  to  0.1  range. 

Figures  6 and  7 are  perspective  drawings  of  the  cockpit  as  seen  from  the  pilot’s  nominal 
viewing  position  and  direction.  Figure  6 shows  entry  points  on  the  lower  front  and  front  side 
surfaces.  Figure  7 shows  primary  reflections  on  the  right-hand  side  of  the  canopy.  Figure  8 shows 
reflection  points  where  the  pilot  has  shifted  his  viewing  direction  20  degrees  to  the  right.  (Note 
that  some  reflection  points  are  shown  on  the  canopy  frame.  This  is  because  the  structure  outline 
is  drawn  as  straight  line  members  between  the  corner  vertices  instead  of  the  curvilinear  members 
used  in  the  computations.  See  Methods.) 


Figures  9 through  13  show  reflection  points  generated  on  one  canopy  surface  by  external 
rays  entering  another  surface.  Figure  9 shows  reflection  points  on  the  top  surface  generated  by 
entry  points  on  the  right  rear  side  surface.  Figure  10  shows  reflections  on  the  top  surface 
generated  by  entry  points  on  the  right  forward  side  surface.  Figure  11  shows  reflections  on  the 
right  rear  side  due  to  entry  points  on  the  left  rear  side.  Figure  12  shows  reflections  on  the  right 
forward  side  due  to  entry  points  on  the  lower  front  surface.  Finally,  Figure  13  shows  reflections 
on  the  right  rear  side  due  to  entry  points  on  the  left  forward  side. 
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Figure  2.  Top,  side  and  front  views  of  canopy  frame. 


Figure  3.  Perspective  view  of  the  cockpit  interior  from  the  pilot’s  nominal  viewing  position  and  direction. 


ray  positions  generating  primary  reflections  on  the 
side  of  the  canopy  as  seen  from  the  pilot’s  position 


Figure  6.  Perspective  view  of  entry  ray  positions  for  the  pilot’s  nominal  viewing  direction  and  position. 


^ Perspective  view  of  primary  reflection  points  for  the  pilot’s  nominal  viewing  direction  and  position. 
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f’igure  8.  Perspective  view  of  primary  refiection  points  for  the  pitot  viewing  20-degrees  to  the  right  side. 
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ire  9.  Entry  ray  positions  on  the  right  rear  side  canopy  surface  and  their 
corresponding  reflection  points  on  the  right  side  of  the  top  surface. 


; 1 0.  Entry  ray  positions  on  the  right  forward  side  canopy  surface  and  their 
corresponding  reflection  points  on  the  right  front  of  the  top  surface. 


ntry  ray  positions  on  the  left  rear  side  canopy  and  their  corresponding 
reflection  points  on  the  top  edge  of  the  right  rear  surface. 


positions  on  the  lower  front  canopy  surface  and  their  corresponding  reflection 
nts  on  the  upper  rear  corner  of  the  right  forward  side  surface. 


[ntry  ray  positions  on  the  left  forward  side  canopy  surface  and  their  corresponding 
reflection  points  on  the  upper  front  edge  of  the  right  rear  side  surface. 


Figures  14  through  19  are  similar  drawings  for  the  copilot’s  position.  Figure  14  shows  “dots” 
for  the  entry  points  of  external  rays  generating  primary  reflections  seen  from  the  copilot’s 
position.  The  figure  shows  that  entry  points  occur  on  the  forward  side  surfaces.  Figure  15  shows 
that  the  corresponding  primary  reflections  occur  on  (1)  the  upper  corner  of  the  lower  front 
surface,  (2)  the  upper  edge  of  the  forward  side  surfaces,  and  (3)  the  upper  front  surface.  The 
associated  reflectance  values  range  In  value  from  0.01  to  0.1.  Figure  16  is  a perspective  drawing 
of  the  cockpit  from  the  copilot’s  position.  The  drawing  shows  reflection  points  where  the  copilot 
has  shifted  his  viewing  direction  45  degrees  to  the  right. 

Figures  17  through  19  show  pairings  between  entry  points  and  reflection  points  by  canopy 
surfaces.  Figure  17  shows  reflections  on  the  upper  front  surface  generated  by  entry  points  on 
the  right  forward  side  surface.  Figure  18  shows  reflections  on  the  lower  front  surface  due  to 
entry  points  on  the  right  forward  side.  Finally,  Figure  19  shows  reflections  on  the  right  forward 
side  due  to  entry  points  on  the  left  forward  side. 


FURTHER  RESEARCH 

The  following  additional  analyses  are  to  be  conducted  for  further  development: 

1.  Investigate  the  reflections  generated  during  realistic  nighttime  lighting  situations  and 
approach.scenarios. 

2.  Investigate  the  solar  glint  generated  as  a function  of  observer-aircraft-sun  angles. 

3.  Investigate  the  use  of  optimization  techniques  to  design  an  optimum  configuration 
for  the  transparent  surfaces  of  the  canopy  to  minimize  both  glint  and  glare  reflections. 


CONCLUSION 

A computer  program  developed  by  HEL  tC'  show  internal  cockpit  reflections  of  external 
point  light  sources  has  been  applied  to  the  Model  YAH-64  Advanced  Attack  Helicopter  (low  glare 
canopy  design).  The  results  show  that  during  nighttime  operations,  ground-light  reflections  are 
possible  on  the  transparent  surfaces  of  the  canopy.  Reflections  are  possible  from  the  top  and  side 
canopy  surfaces  for  the  pilot  and  the  front  and  forward  side  surfaces  for  the  copilot.  The  results 
are  an  improvement  over  the  flat  plate  canopy  design  since  reflections  are  limited  to  certain 
portions  of  these  surfaces.  Where  reflections  actually  occur  depend  upon  the  particular  lighting 
situation  and  flight  scenario. 
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ry  ray  positions  generating  reflections  on  the  right  hand  side 
canopy  as  seen  from  the  copilot’s  nominal  position. 


Figure  16.  Perspective  view  of  primary  reflection  points  for  the  copilot’s 
viewing  45-degrees  to  the  right  side. 
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Figure  1 7,  Entry  ray  positions  on  the  right  forward  side  canopy  surface  and  their  corresponding 


i.  Entry  ray  positions  on  the  right  forward  side  surface  and  their  corresponding 
reflection  points  on  the  upper  right  corner  of  the  lower  front  surface. 


Figure  19.  Entry  ray  positions  on  the  left  forward  side  surface  and  their 
corresponding  reflection  points  on  the  upper  edge  of  the  right  forward  side. 
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APPENDIX  A 

INTERNAL  REFLECTIONS  FROM  CYLINDRICAL  SURFACES 


The  ray-tracing  techniques  used  in  the  computation  are  those  of  Appendix  A,  reference  5 
with  appropriate  augmentation  for  computing  the  intersection  of  a straight-line  ray  with  a 
cylindrical  surface.  As  in  reference  5,  coordinates  are  measured  in  a cartesian  coordinate  system 
(x,y,z)  with  the  y-axis  along  the  longitudinal  axis  of  the  aircraft,  the  z-axis  directed  into  the 
canopy  top  and  the  x-axis  orthogonal  to  the  other  two  axes. 

A straight-line  ray  is  specified  by  an  origin  point  and  its  directional  cosines 

(aj.bj.Cg).  We  wish  to  compute  the  intersection  point  of  this  ray  with  a cylinder  and  determine  if 
the  intersection  point  is  enclosed  by  the  edges  of  a surface  on  the  cylinder.  A cylinder  is  specified 
by  three  parameters  (1)  an  origin  point,  Pq  (xQyQ,ZQ)  on  the  cylindrical  axis,  (2)  directional 
cosines  (aQ,bQ,CQ)  for  the  axis,  and  (3)  its  radius,  r^..  A cylindrical  surface,  then,  is  specified  by 
these  cylindrical  parameters  and  a consecutive  sequence  of  N adjacent  corner  vertices  (Pj 
(xj,yj,Zj),  I = 1,N)  on  the  cylinder  and  their  coordinates  in  the  rectangular  space. 


(1 ) Intersection  of  a straight-line  ray  with  a cylinder 


The  intersection  point  of  a line  with  the  surface  of  a cylinder  may  be  computed  given 
the  origin  (x5,y5,Zg)  and  directional  cosines  (a5,b5,Cg)  of  the  line,  and  the  origin  (xQ,yQ,Z(j)  and 
directional  cosines  (aQ,bQ,CQ)  of  the  cylinder’s  axis  and  the  radius,  r^.,  of  the  cylinder.  Note  that 
the  coordinates  of  the  intersection  point,  Pj  (x|,y|,Zj) , are  given  by 

Xj  = Xj  + aj  • Rj, 


yi  = ys+bs-Rs> 

^i  ~ ^s  ^ S *^s' 

where  R^  is  the  length  along  the  line  between  its  origin  and  the  intersection  point. 


0) 


Consider  now  another  line  from  the  intersection  point,  Pj  orthrogonal  to  the 
cylindrical  axis.  The  dot  product  of  the  unit  vectors  must  be  equal  to  zero;  i.e.. 


ao  • (Xj-xJ  + bg  • (yj-y^.)+  c^  • (zj-z^.)  = 0,  (2) 

where  the  point,  P^.  is  at  the  intersection  of  the  axis  and  this  line.  Note  that  the 

coordinates  of  this  point  are  given  by 
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(3) 


~ ^0  ^ ^0^0' 


yc^Vo’^^o'^o’ 


~ ^ *'0^0’ 


where  Rq  is  the  distance  along  the  cylindrical  axis  between  the  points  Pq  and  P^,. 

Substituting  equations  (1)  and  (3)  into  equation  (2),  we  obtain  a relationship  between 
the  two  distances, 


Ro  = a'^s'^3, 


(4) 


where 


6 = ao(Xs-Xo)  + b^iy^-yj  + Co(z5-Zo). 


Note  that  the  parameter  a equals  the  dot  product  of  the  unit  vectors  for  the  straight-line  ray  and 
the  cylindrical  axis.  The  parameter  6 equals  the  length  of  the  projection  of  the  line  P^P^  onto  the 
cylindrical  axis. 

The  line  connects  a point  on  the  surface  of  the  cylinder  and  one  on  the  axis.  Since 
it  is  orthogonal  to  the  axis,  its  length  equals  the  radius  of  the  cylinder  and 

2_ 


rc  =(xi-Xc)^+(yj-yc)^+(Zi-Zc)' 


(5) 


Substituting  equations  (1),  (3)  and  (4)  into  equation  (5)  produces  the  quadratic  expression 


a*  + 2bRj  + c-r^.^  = 0, 


(6) 


where  a = a-|  ^ + 32^  + 33^, 

b = a^b-|  + 32b2  + a3b3, 


c = bi^  + bo^  + b 


3 ' 
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:emm 


and  a-|  = «aQ-  a^, 


bi  =Xo-Xs  + 0ao, 


a2  = «bo-b5, 

b2  = yo-  Vs  + ^bo- 
^3  "“‘'0  “ S' 


h = ^o"h'^^^o- 


Solving  equation  (6)  leads  to 


(7) 


where  the  sign  of  the  square  root  is  determined  by  the  restriction  that  Rg  be  larger  or  equal  to 
zero,  Rg>  0.  Using  equation  (7)  in  equation  (1)  leads  to  the  coordinates  of  the  intersection  point 

Pi- 


(2)  Determine  if  the  intersection  point  is  enclosed  by  the  edges  of  a cylindrical  surface. 

The  techniques  of  section  (5),  Appendix  A,  reference  5,  can  be  used  to  determine  if 
the  intersection  point  of  a ray  is  within  a convex  planar  figure  with  straight  edges.  A 
transformation  of  the  surface  of  the  cylinder  into  a plane  will  allow  the  application  of  these 
techniques.  We  assume  that  the  curvilinear  edges  are  so  designed  that  they  transform  into  straight 
lines  (see  Reference  1). 

An  appropriate  transformation  is  one  that  will  change  the  cartesian  coordinates  of  a 
point  into  coordinates  of  a cylindrical  coordinate  system  A suitable  cylindrical  system 

could  be  defined  as  follov/s.  The  ? coordinate  is  the  distance  along  the  cylindrical  axis  measured 
from  the  origin  point,  Pq.  Positive  values  would  be  in  the  direction  of  the  axial  unit  vector.  The  C 
coordinate  is  the  distance  along  the  cylindrical  surface  measured  in  a plane  normal  to  the 
cylindrical  axis.  The  distance  is  measured  from  a reference  plane  containing  the  cylindrical  axis 
and  a predefined  reference  vector.  Positive  values  would  be  counterclockwise  displacements  as 
seen  facing  the  axial  unit  vector  for  a right-handed  coordinate  system.  The  n coordinate  is  the 
distance  above  or  below  the  cylindrical  surface  measured  along  a normal  to  the  surface.  Positive 
values  are  above  the  surface  away  from  the  cylindrical  axis. 
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We  may  compute  the  cylindrical  coordinates  ^ ,C;i)  of  a point  P from  the  rectangular 
coordinates  (x,y,z)  by  first  constructing  a line  from  the  point  P orthogonal  to  the  cylindrical  axis. 
The  rectangular  coordinates  (x  cjYc^c^  **c  intersection  of  the  two  lines  are 

given  by  equation  (3)  in  section  (1)  above.  The  dot  product  of  the  unit  vectors  for  the  two  lines 
is  equal  to  zero  and  the  distance  Rq  along  the  axis  between  the  origin,  Pq,  and  the  corresponding 
intersection  point  is: 


Rq  = ao(x-Xo)  + bo(y-yo)  + Cq(z-z^).  (8) 

The  distance,  R,  of  the  line  P^-P  is  given  by 

R = i(x-xj2  + (y-y^)2  + (z-zj2('/2  (9) 

where  equation  (8)  is  substituted  in  equation  (3)  for  the  coordinates  The  directional 

cosines  of  the  line  are  given  by: 


ac  = (x-X(.)/R, 


be  = (y-yc)/R. 


(10) 


Cc=  (z-Zc)/R, 

We  next  construct  a reference  vector  for  the  cylindrical  coordinate  system  as  follows. 
Let  the  vector  be  normal  to  the  cylindrical  axis  and  in  a vertical  direction.  The  z-coerdinate  of 
the  point  Pq  is  greater  than  zero,  Zq>  0,  and  we  consider  the  truncated  point,  Pq*  (x^  ,yQ,0),from 
which  a line  is  constructed  orthogonal  to  the  cylindrical  axis.  Let  the  line  and  axis  intercept  at 
the  point  Pjj  (xj,,yjj,Zjj).  Then  the  distance  along  the  axis  from  Pq  to  Pj,  is  R|j  = - CqZq  as  given 
by  equation  (8).  The  coordinates  of  point  Pjj  are 

^b  ~ ^o~^o''o^o> 


Yb  ^0  bo^o^o’ 


Zb  = (1-Oz 


o> 


as  given  by  equation  (3).  Also,  the  length  of  the  line  Pq  P(j  is  R 


= zjl- 


from  equation  (9). 


Use  was  made  of  the  fact  that  the  sum  of  the  squares  of  a set  of  directional  cosipes  is  equal  to 
unity;  i.e.,  aQ^  + bQ^  + Cq^  = 1.  And  finally,  the  directional  cosines  of  the  line  pjpjj  are  given  by 
equation  (10)  as 


^R  ^0*^0 


bR'-bpCp, 

Vl-Cp^ 


(11) 
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These  are  the  directional  cosines  of  the  reference  vector  contained  in  the  reference  plane  for  the 
cylindrical  coordinate,  ?. 

We  construct  a vector  with  directional  cosines  (ap,bp,Cp)  normal  to  the  reference  plane 
from  the  cross  product  of  the  axial  unit  vector  and  the  reference  vector;  i.e., 


^p'^o'^R  ^R^'o’ 
bp  = -(aoCR-aRCo). 


Cp  = »o^R~^R^o- 

This  vector  at  the  cylindrical  surface  Is  in  the  direction  of  decreasing  values  of  the  C-coordinate. 

We  let  the  angle  '/»q  be  the  angular  displacement  of  the  line  P^P  from  the  reference 
plane  The  line  pTp  is  contructed  from  the  point  P and  is  orthogonal  to  the  cylindrical  axis.  The 
angle  I'o  is  measured  in  a plane  normal  to  the  axis.  It  is  equal  to  the  arccosine  of  the  dot  product 
of  the  unit  vextors  of  the  line  p?  and  the  reference  vector;  i.e. 

(13) 

ijiQ  = arccosine  (3r  * a^.  + bj^  • b^.  + cj^  • c^) 

Note  that  the  dot  product  of  the  unit  vectors  for  the  line  P^P  and  the  reference  plane  normal, 
i.e., 

(14) 

Q = 3c  • I’c  ’ ^p  ^ ^c  ' ‘'p 

has  the  same  sign  as  does  the  angle.  That  is,  if  Q < 0,  then  < 0,  and  if  Q > 0,  then  % > 0. 


The  cylindrical  coordinates  of  the  point  P are  determined  as  follows.  The  ^-coordinate 
is  the  distance  along  the  cylindrical  surface,  measured  in  a plane  normal  to  the  cylindrical  axis, 
from  the  reference  plane  to  the  line  F^.  The  distance  is  given  by 


where  is  the  angular  separation  of  the  P^-P  given  by  equations  (13)  and  (14).  The? -coordinate 
is  the  distance  along  the  axis  from  the  origin,  Pq,  to  the  point,  P^.,  i.e., 

(16) 


where  R^  is  given  by  equation  (8).  Finally,  th_en -coordinate  is  the  distance  of  the  point  P above 
or  below  the  cylindrical  surface  along  the  line  i.e.. 


n =R-r, 


(17) 


where  R is  given  by  equation  (9). 
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We  may  use  the  techniques  of  section  (5),  Appendix  A,  reference  5 by  first  applying 
the  transformations  of  equations  (15),  (16)  and  (17)  to  be  points  of  interest.  These  are  (1)  the 
ray  origin,  Pj,  (2)  the  intersection  point,  Pj  of  the  ray  with  the  cylinder,  and  (3)  the  corner 
vertices  of  the  cylindrical  surface,  (Pj,j=l,N). 


(3)  Computing  transmitted  and  reflected  rays 

The  reflected  and  transmitted  ray  components  and  their  associated  values  of 
transmittance  and  reflectance  can  be  computed  using  the  techniques  of  section  (6),  Appendix  A, 
reference  5.  The  transparent  surfaces  of  the  canopy  are  assumed  to  have  infinitesimal  thickness. 
All  internally  refracted  rays  leave  at  the  same  surface  point  as  does  the  initially  transmitted  ray, 
and  this  point  is  the  same  as  that  at  which  the  incident  ray  reaches  the  surface  (see  references  2 
and  3.) 


The  angle  of  incidence,  0q,  between  the  incident  ray  and  the  surface  normal  is 
determined  by  the  dot  product  of  the  unit  vectors;  i.e., 

0Q  = arccosine  (aj  • a^  + bj  • b^^  + Cg-Cj,),  (18) 

where  ag,bg,Cg  are  the  directional  cosines  of  the  incident  ray,  and  ap,,bp,,Cp  are  those  of  the 
surface  normal. 

The  normal  to  the  cylindrical  surface  is  measured  at  the  intersection  point,  Pj,  of  the 
inddent  ray.  The  normal  is  directed  into  the  cockpit  volume  by  convention.  It  is  along  the  line 
of  section  (1 ),  and  its  directional  cosines  are  given  by: 

^n  ~ 

^n  = (Vc-Vi)/ '"c'  (^9) 

Cn  = (Zc-Zj)/rc. 

Note  that  the  coordinates  Xj.,yj.,Zj,  of  the  point  P^.  are  given  by  equation  (3)  while  those  Xj,yj,Zj 
of  Pj  are  given  by  equation  (1). 
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APPENDIX  B 


COMPUTER  PROGRAM 


The  computer  program  is  programmed  for  a disk-based  batch-operating  system  in 
FORTRAN  IV  language  (see  references  4 and  7).  The  program  is  as  listed  in  Appendix  C, 
reference  5,  except  for  a few  additions.  A short  list  of  additional  subroutines  follows.  The 
program  is  attached. 

1.  CALC— Controls  the  computation  of  the  reflection  point  for  an  incident  ray  and  stores 
the  calculations  on  file,  if  any.  Called  by  CONTL. 

2.  INTCY— Computes  the  intersection  point  of  an  incident  ray  with  a cylinder  and  tests  the 
point  for  enclosure  within  a cylindrical  surface.  Called  by  CALC. 

3.  TRSCY— Transforms  rectangular  coordinates  into  cylindrical  coordinates.  Called  by 
INTCY. 

4.  PREC— Draws  perspective  of  cockpit  as  a function  of  viewing  angles,  and  shows  entry  and 
primary  reflection  points.  Called  alone. 

5.  DRWCT-Draws  distribution  of  points  on  all  surfaces.  Called  alone. 

DRWCP-Draws  side,  top  and  front  views  of  canopy  and  entry  and  primary  reflection 
points.  Draws  pairing  of  entry  and  reflection  points  by  surface.  Called  alone. 

6.  DRWCF-Draws  sides,  top  and  front  views  of  canopy.  Called  by  DRWCT  and  DRWCP. 
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/NEkJCG(HCAAH) 

/JCe(HCAAH) 

/CREATE(FINKAH> 

/FORT 

C MAINLINE,  RA>  TRACING  FOR  PILCT 

C 3 OIHEKSIONAL  FLAT/CYLINORICAL  CANOPY  NITH  OBSTRUCTIONS 
SUBRCLTINE  CCNTL 
CCMMCN/FCRG/LFG 
COMMON/GAREA/IDFIL(6COO) 

COMHON/ANG/ANN.BNN 
CATA  ANN,BNN/S5.,S0./ 

CALL  GETDEV(LLN,'FINKAH* ,10) 

CALL  GETDEV(LFG,'FORGRA*  ,-l) 

CALL  CAND 
CALL  CANLILLM 
CALL  CANT 

CALL  OEVTaLN,C,U) 

RETLRN 

ENC 

/FORT 


CONTROLS  INPLT  OF  DATA 
SLCRCLTINE  CANO 
DIMENSION  AN(2) 

CATA  AN/2FYS,2HN0/ 
hRITEdtSSG) 

•998  F0RMAT(2X, 'CANOPY  DATA  ON  FILES  (YS  OR  NOP) 

REACH, ICCDA 
IF(A.EQ.AMl))  GO  TO  1 
CALL  REAOV 
CALL  NORKL 
GO  TO  2 

1 CONTINUE 
CALL  REAOF 

2 CONTINUE 
MRITEtUlCCC) 

1000  F0RHAT(2Xf 'PRINTOUT  CANOPY  DATA!  (YS  OR  NO)') 
RGAD(1,1CC1)A 
ICCl  F0RNAT(1A2) 

IF(A.EQ.AN(1) ) CALL  TELTY 
kRITE(l.UC2) 

1002  F0RMAT(2X, 'DISPLAY  CANOPY  CONFIGURATIONS  (YS  OR  NO)') 
REACH, ICCDA 

IF(A.EQ.AN(2))  GO  TO  3 
CALL  DRkCN 
kRITEH,lC03> 

1003  F0RMAT(2X,'HARC  CCPYS  (YS  OR  NO)') 

REACH, ICCDA 

IF(A.EQ.AN(2) ) GO  TO  3 
hRlTE(l,lCOA) 

lOCA  F0RHAT(2X, 'PRINTER,  TURN  ON') 

PAUSE 

CALL  GS3i(l) 

PAUSE 

3 CONTINUE 
CALL  GHLT 
)|RITE(1,1CCS) 

IOCS  F0RHATI2X, 'DISPLAY  PERPECTIVES  (YS  OR  NO)') 

REACH, ICCDA 
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IF(A.EQ.4N(2) ) GO  TC  4 
CALL  PERP(C) 
bRITE(ltlCC3) 
REAOdtlCCDA 
IF(A.EQ.AR(2nGGTG  A 
CALL  GS31(1) 

PALSE 
A CONTINUE 
CALL  GHLT 
RETURN 
ENC 

/FORT 


jffi 


CONTROLS  CALCLLATION  CF  REFLECTION  POINTS 
sueRCuriNE  canulln) 

COPPCN/CAN/NT*N8,NC,KP,NA,NO.NV(100)»PXV(lOO,8),FYV(100>8)tPZVaOO 

Qt8) 

CIHENSICN  AN(2) 

CATA  AN/2HYS.2HNO/ 

CATA  AIl,AI2/2.,2./ 

CALL  SEEK(LLN«C) 

REACILLNtNVN 

NRITE(l.SSS) 

999  FaRP.AT(2X,'CCHPLT£  PRIHARY  REFLECTION  POINTS!  (YS  OR  NO)*) 
REACdtlCCDA 
lOGl  FCRPATdA2) 

IF(A.EO.AN(2))  GO  TG  300 
NVLL»C 

CALL  SEFK(LLNfNVh) 

NRITE(LLN)NYLL 

INDEX  PILCT  VIEWING  DIRECTION 

Al  ELEVATION,  ANGLE  PROP  2 AXIS  TOWARD  X AXIS  IN  X 2 PLANE 
A2  A2IHUTH,  ANGLE  FROP  X 2 PLANE  TOWARD  Y AXIS  . 

2 AXIS  TOWARD  UPWARD,  X AXIS  TOWARD  LEFT  FACING  FRONT  OF  CANOPY  AND  Y AXIS 
TOWARD  RACK  OF  CANOPY  ALONG  LONGITtCINAL  AXIS  OF  AIRCRAFT 
A1«C. 

10  A2«AI2 

11  A2»A2-AI2 
IF(A2.GE.>9G.)  GO  TO  15 
Al>Al«Ad 

IFIAl.GT.ieC.)  GO  TC  299 
WRITEd,99!)Al 
99S  F0RNAT(2X,F1C.A) 

GO  TO  10 
IS  CONTINUE 
IS«0 

CALL  CALC(A1,A2,IS,LUN,NVLL) 

IFdS.EQ.NAICALL  CALCIAl  ,A2,NA,LUN,NVLL) 

CO  TC  ll 

299  CONTINUE 

CALL  SEEK(LLN,NVW) 

WRITE(LLN)NVLL 

300  CONTINUE 
RETURN 
ENC 


CONTROLS  OUTPUT  CF  CALCULATIONS 


35 


o n o 


SU6RCLTINE  CANT 
CINENSICN  AN(2) 

CATA  AN/2HYS,2HN0/ 

(|RITE(1«1CC6) 

1C06  F0RNAT(2X,*PERPEC7IVE*/2X,*BAY  ENTRENCE  P0IMS*/2X, 'RIGHT  HAND  SIE 
CE  COCKPIT  IY£  CR  NO)') 

REACdtlCCl  )A 

lOCl  FORMATI1A2) 

IFIA.EQ.ANC2) ) GO  TO  40 
CALL  PERP(l) 
hRlTE(UlCC3) 

1CC3  F0RHATJ2X, 'HARD  CCPY+  <YS  OR  NO)') 

REAOd.lCCDA 
IF(A.Eg.AN(2) )GCTC  40 
CALL  GS3K1) 

PALSE 

40  CONTINUE 
CALL  GHLT 
WRITEd.lCC?) 

ICO?  F0BMAT(2X,'P£RPECTIVE'/2X,'PR1HARY  REFLECTIONS!') 

REAOdtlCCDA 
IF(A.EQ.AN(2))  GO  TO  41 
CALL  PERP(2) 

WRITEd«lCC3) 

REAOdtlCCDA 
IF(A.EQ.AN(2) )GOTC  41 
CALL  GS3K1) 

PALSE 

41  CONTINUE 
CALL  GHLT 
RETLRN 
ENC 

/FORT 


PRIMARY  REFLECTION  ONLY,  SURFACE  NK  TRANSMITTER  ONLY 
SU8RCLTINE  CALC  (A1  »A2  *NK  ,LLN  »NVLU 

COMMON/C  AN/NT, NBtNCtNPtNA, NO  »NVdOO ) tPXVUOOy  8 ),PYVdOO»8)»PZVdCO 

Qf8) 

C0MMCN/PILC7/XC»YCfZ0 
COMMON/LINE/AS«BS*CS»XS*YSyZStAC»BC»CC 
CAT!  PI/3,14159/ 

AlM»Al«PI/ieC. 

A2M-A2«PI/18C. 

AS«C0SfA2P)*SIN(AlM) 

BS»SIN(A2M) 

CS»C0S(A2M)*CCS(A1M) 

R«l. 

XS«XO 

YS*YO 

ZS»ZC 

IPLN*C 

INK=C 

16  CONTINUE 

CO  20  1S=1,NC 
IF(IS.EQ.INK)GC70  20 
ISK  = C 

IF(IS.LE.NA)CALL  INTEC dSK, IS, XR,YR,ZR) 

IF( IS.G7.NA)CALL  I NTC Y C I SK , I S , XR , YR,ZR) 

IF(  ISK.GT.OGOTO  25 
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2C  CONTINUE 

return 

25  CONTINUE 

IFUS.LE.NORETURN 
CAUL  COMP(ANG.RT,TT) 
T-TT^R 

IFUS.EQ.NK  >GCT0  AO 

IPUN«IPUN4l 

AI»-AS 

BI»-BS 

CI»~CS 

RX»AS*AC4eS*BC4CS*CC 
AR»-AS42«ARX*AC 
BR»-BS42.*RX*9C 
CR«-CS42.*RX*CC 
1F(IPUN.EQ.2»GCT0  3C 
IFIR.LT..CCCCURETURN 
IPP*IS 
ASP«AS 
BSP’BS 
CSP»CS 
XP«nR 
YP»YR 
2P«ZR 
AIP>A1 
elP*Bl 
C1P»CI 
ARP«AR 
BRP«6R 
CRP»CR 
AS—AR 
es«-BR 
CS*-CR 
XS*XR 
YS»YR 
2S»ZR 
6C  TO  18 
30  CONTINUE 

IF(T*l.T..CCCCl)REURN 

Q/l) 

NK«IPP 
RETURN 
AC  CONTINUE 
INX’NK 
GO  TO  18 
ENC 


P.FAC  IN  SURFACE  VERTICES 

S2!l5r«UM“B.KC,KP.N..KC,xv(ioo).p»v.ioo,»>.PYv<ioo.e..Piy<ioo 

OfO 
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COMMON/VERT/)(V(2CC)«YV(2CC}  ,ZV(200},^VR(100f8) 
C0RHQN/CYL/NCY|NSC(10)  ,NSP(10,10),XC(10),YC(10)»ZCaO)»«E(10),eE(l 
QC)tCEflO)*RC(lC) 

COPMON/PILCT/XC*YC>ZO 

REA0<2tlCCC) 

ICCO  FORHATf) 

KEAO(2*lCCl)NTtNBtNCtNP*NA»hD 
ICCl  FCRMATI/6(2X*I3)) 

READ(2«1CCC) 

REAO(2flCC2)(NV(n«I>t,NOI 
1CC2  F0RPAT(16(2X,I3) I 
REA0(2«1CCC) 

CO  10  J«lt8 

REA0(2»l€C2)(NVR(ItJ)«I«l«KD) 

10  CONTINUE 
REAC(2«1CCC) 

REA0(2tlCC3)(XV(It«I>ltNTI 
1CC3  F0RNAT(8(2X«F7.A) ) 

REA0(2tlCCC) 

REA0(2flCC3MYV(ntI>ltNT) 

READ(2flCCC) 

REA0(2flCC3l(ZV(n*I>l»NT) 

REA0(2«1CCC) 

REA0(2flCC2)NCY 

REA0(2tlCC2MNSC(ntI*lfNCY) 

CO  20  1«1«KCY 
KP«NSCm 

20  REA0(2tlCC2)(NSPntR)tK«l«KP) 

READ(2flC0A)(XCm*YC(n  t ZC  (I ) * AE  ( I)  »BEI  I)  »CE  ( I)  *RC  ( I )•  I-l.NCV ) 
1004  F0RNAT(7(2X«F7,4}  I 
REA0(2tlCCC) 

REAC<2tlCC3)XOtYO«ZO 

RETURN 

END 

/FORI 

C 

C 

C 

C ESTA8USH  SURFACE  NORPAU  FOR  EACH  PLATE  SURFACE 
C SURFACE  NORMALS  DIRECTED  TOWARD  COCKPIT  INTERIOR 
SUBROUTINE  NCRRL 

CONMON/CAN/NT«NBfNCfNP*NA»ND»KV(100),PXVI100t8)»PYV(100,8),PZV(100 

QfBI 

C0NM0N/VERT/XV(20C)tYV(200)«ZVf 200) «AVRf 100*8) 

CQNMON/NORM/AXM  ICO)  tAYN(lCO)  »AZNf  I'JO) 

C0MHQN/CYL/NCYfNSC(10)»NSP(l0«10)»XCa0),YC(10)*ZC(10),AE(10)»BE(l 

00)tCEUC)»RC(lC) 

CONMON/PlLOT/XCtYCtZO 

NVIi*C 

CALL  GETDEV(LUNt*FINKAH«*10) 

CALL  SEEKILUNtC) 

NRlTEaUNINVk 

WRITE (LUN)NT,NB»NC»NP,KA*KO 

CO  10  i>:*N0 

NK*NVm 

WRITEILUN)NK 

CO  5 K*ltNK 

KV«NVR(I*K) 

PXV(1*K)»XV(KV) 

PYV(1*K)-YV(KV) 
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PZV(I»K)=ZV(KV) 

V*RITEtHJN)PXV(I,K)  ,PYVn  ,K)  ,PZV(I.K» 

5 CONTINUE 
K*2 

7 A1=PXV(I,K)-PXV(I,1) 

A2»PXVa,K  + l)-PXVn,l) 

Bl*PYVn,K)-PVV(I  ,1) 
e2»PYV(I,K*l)-PYV(I,l) 

Cl»PZV(I,K|-PZV(l,l) 

C2»PZV(I,K-m-PZV(I,l) 

Pl»SQRTIAl**2*Bl**2*Cl**2l 

A»(Al*A2*ei*B2*Cl*C2)/(Pl*P2) 

IFIABSIAt.LT.l.)  GO  TO  9 
K»K  + 1 

IF(K.EQ.NK)  GO  TO  10 
60  TO  7 
9 AN«ACCS(A) 

R'l./IPIOPZ^SINIANI) 

AXNm*><ei*C2-Cl*B2)*R 
AYN(n=-(Al>»C2-Cl*A2)*R 
AZMU»+(A1*82-A2*B1)*R 
hRITEaUNIAXMIltAYNm  tAZM!) 

10  CONTINUE 

NRITE(LUN)NCY*(NSC(n«I«l*NCY) 

CO  20  I«ltNCY 
KP*NSCn  I 

20  URITE(LUN)(NSP(ItK)*K^l*KP) 

hRnEiLUNHxc(ntYcmtZcmyAEmfBEmfCEm»Rcmti«itNCY) 

NRlTEILlNIXCfYCtZC 

NVN>NEXREC(UN) 

CALL  SEEK(LLNtC) 

NRITEILININVN 
CALL  OEVTILLNfGtO) 

RETURN 

ENC 

/FORT 

C 

C 

C READ  FILE  FOR  CANOPY  DATA 
SUBROUTINE  REACF 

COf<MON/CAN/NT  «NB*NCtNPfNA«NOtNV  1100)  tPXVIlOO  *8 ) »FYV(  100*8  ItPZVI  100 
Q*8) 

COPNON/NCRP/AXNUCO)  .AYN(ICO)  yAZMlOOI 

CONNON/CYL/NCY*NSC(10)*NSPaO«10)*XC(10)*YC(10)  »ZC(10l*AE(10l*eE(l 
QC)«CE(1C)«RC(1C) 

CONHON/PlLCT/XC,YCtZO 
CALL  66T0EVtLUN,»FINKAH»,10) 

CALL  SEEK(LLNfO) 

REAOlLUNINVb 

REAO(LLN)NT»NB«NC»NP*NA*ND 

00  10  I«1*N0 

READILUNINK 

NVin*NK 

DO  S K>1*NK 

READILUNIPXVd  *K)  ,PYV  1 1 »K)  »PZV  (I  *K ) 

5 CONTINUE 

READILUNIAXMI  )»AYN  n.AZMI) 

10  CONTINUE 

READILUNINCY*  (NSCm  tl>l*NCY) 
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DC  2C  I=1«NCY 
KP«NSCm 

2C  REACIUJNHNSP(ItK),K«l*KP) 

REAO(LUNI(XC(I)»YC(I)»ZC(I)fAE(n  teEm«CE(ntPC(I)«I>ltNCY) 

REAC(LUN)XOfYOtZO 

CALL  DEVT(LlNtCtO) 

RETLRN 

ENC 

/FORT 


PRINTOUT  OF  CANOPY  DATA  FOR  REVIEh 
SUBROUTINE  T6LTY 

CORHON/CAN/NT«NBtNCfNPtNA»KD»NV(lOOI,PXV(100t8)  tPYV(100»8l»PZVUOO 
Q»8) 

C0NH0N/VERT/XV(20C) »YV(200) tZV(200)tNVR(100f8) 
C0HPCN/N0RP/AXN(1C0)«AYN(1C0) »AZN(100I 

C0PM0N/CYL/NCY«NSC(10) *NSP (10 tlO ) »XC (10) , YC (10 ) t ZC I 10 )» AEIIO }» BE( 1 
QC)»CE(IC),RC(1C) 

COPMON/PILCT/XC,YC,ZO 

RRITEdfSSB) 

S98  F0RHAT(2Xt'PRINTERt  TURN  CK*) 

PAUSE 

)iRITE(3.IC0C)NT,NBtNC»NP*KA,KC 
loco  F0RMAT(18(2X,I3)) 

99«  FORHATO 

WRITE(3f997) 

997  F0RMAT(/2X,»VERT£X  POSITION  DATA') 

NRITE(3«lCCU(XV(I)tI>l*KT) 

NRITE(3t996) 

kRITE(3fiC01)  (YVm»I^ltNT) 
kRITE(3f9S6) 

kRITEOtICOl)  (ZV(I)«I*1»NT) 
kRITE(3f996) 

kRITE(3f 1CCC)((NVR(I tJ)»I«lfNC)tJ«l»8) 
kRITE(3»996) 

CO  10  I*l«NO 
KN*NV(I ) 

kRITE(3t995)( fKN 

995  F0RMAT(/2X»'SU'RFACE*fI3,*N0.  VERTICES*  1 13) 

CO  5 K*lfXN 

kRlTE(3,lCCl)PXV(I,K)*PYV(I,K),PZV(I,K) 

1001  F0RHAT(2Xt3(F7.2»2X) ) 

5 CONTINUE 

kRITE(3,lC02)AXN(I)«AYN(I)*AZMI) 

1002  F0RMAT(7(2X,F1C.9)) 

1C  CONTINUE 

kRITE(3«99A) 

99A  F0RPAT(/2X, 'CYLINDRICAL  DATA*) 
kRITE(3»lC0C)NCY 
kRITE(3,lCCC)(NSC(I ),I»1,NCY) 

CO  20  I’ltNCY 
KP«NSC(I) 

20  kRITE(3.1CCC)(NSP(IfK).K»l,KP) 

kRITE(3»lCC2)(XC(I)«YC(I) »ZC(I),AE(I)»BE(I),CE(I)»RC(I)»I^1»NCY) 
kRlTE(3«993) 

993  F0RPAT(/2X, 'PILCT  POSITION*) 
kRITE(3tlCC2)X0tYC»Z0 
RETURN 
ENC 
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/FORT 


DETERMINES  IF  RAY  STRIKES  CONVEX  SURFACE 
SUBRCtTINE  INTEC( ISK,IS»XR,YR*ZR| 

CCMMCN/CAN/NT.NB,NCtNP*NAtND»NV(100)tPXV(100*8)*PYV(100,8}»PZVnOO 

Q*8) 

COMMON/NCRM/AXMICO)  tAYNIlCO)  tAZN(lOO) 
COMMON/LINE/AS«BS»CS«XStYS*ZS»AC»BC,CC 
CK*AXN(IS)*AS^AYN(1S)«BS4AZMIS)«CS 
IFICK.GE.C.)  RETURN 

RAY  STRIKES  SURFACE  IN  OLTKARC  DIRECTION 
I»1 

IF(XS.EQ.PXV(IS«I  ).ANO.YS.EC.PYV(IS»n.ANC.ZS.EC.PZV(IStini«I«l 
S*(AXN(IS)*(PXV(IS»I)-XS)«AYN(IS)*IPYVIIS«n-YS)«AZMISi*(PZV(IS*I 
c)-zsn/CK 
XR«AS«S*XS 
YR«BS>*S+YS 
ZR«CS*S4ZS 
A1«PXV(IS,1)-XS 
Bl»PYV(IStn-YS 

ci«pzv(is»n-zs 

Pl»XR-XS 

P2«YR-YS 

P3*ZR-ZS 

1N>NV(1S) 

CO  IC  I«1*IN 
IC»l4l 

IF(I.Eti.IN)  IC«1 
A2*PXV(ISiIC)«XS 
B2«PYVnSfIC)-YS 
C2«PZV(IS,IC»-ZS 

Q»Pl*IBl*C2-B2*Cl»-P2*(Al*C2-CMA2)*P3*(Al*e2-Bl*A2I 

IF(Q.GT.C«)RETIRN 

RAY  STRIKES  SURFACE  ON  ENCLOSED  SICE  OF  SURFACE  ECCE 
Al»A2 

el«B2 

Cl*C2 

10  CONTINUE 

RAY  STRIKES  ENCUCSED  SURFACE 
ISK^l 

AC«AXN(IS) 

8C*AYN(IS) 

CC*AZN(IS) 

RETURN 

END 


CCMPUTES  INTERSECTION  POINT  OF  LINE  WITH  CYLINDER 
SUBROUTINE  INTCY ( I SK t IS t XR » YR »ZR ) 

COMMON/CAN/NT  *NB»NCfNP»NA«NCtNV  (100  )tPXVa00»8)»PYV(10O>8)fPZV(  100 
QtB) 

CCMMON/LINE/ASfBStCSf XStYSfZStACtBCtCC 

COHMON/CYU/NCYtNSC(lO) »NSP ( 10 *10 ) »XC (10 ) *YC ( 10 ) t ZC( 10 ) t AE( 10  It  BE ( 1 
CC)«CE(LCI»RC(1C) 

CCMMCN/PILCT/XP^YP,ZP 
C DETERMINE  CYLINDERICAL  SURFACE  WHICH  CONVEX  SURFACE  IS  A PART  OF 
CO  2 1 = 1, NCV 
KP=NSC(I » 
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CO  2 K=l,KP 

IF(NSP(I,K).EQ.IS»GCTO  5 
2 CONTINUE 
RETURN 
5 CONTINUE 

C DETERMINE  INTERSECTION  POINT 
XC*XC(I) 

YC«YC(n 
ZO>ZC(I ) 

RC»RC(I) 

AO>AE<I) 

BO»BEm 

co*cEm 

ROS»SQRT(  (XS-XC)»*2-MYS-Y0)A»2*(ZS-Z0)*A2» 

ACC*IXS-XC)/ROS 

BCC»(YS-YC)/ROS 

CCC»(ZS-ZC»/RCS 

A>AO*AS4BC«BS-*>CO«CS 

B»RCS*(A0*ACC4BC*80OKO*CCO) 

Al*AS-A*AC 

A2»RS-A*BC 

A3»CS-A*CC 

Bl*ACC*RCS>B’»AO 

B2«B0C*R0S-e*BC 

B3«C00*RCS-B*CC 

A«Al**2+A2**2*A3**2 

B«Al*Bl4A2*fl2*A3*B3 

C»B1A*2*B2<"*2*B3**2 

ABa-B/A 

B0a(  SQRT(B**2-AA{C-R0**2)  n/A 
RS*AB>BB 

IFIRS.LT.C. )RSaAB«BB 
XR*XS4A$«RS 
YR»YS*RS*RS 
ZR«ZS+CS*RS 

R«RS*(AO*AS*eO*BS+CC'»CS»  + POS*(AO*AOO^BOABOO»CO*CCO  » 

X1«X0«A0AR 

Yl»YC*RCAR 

Zl«ZO«CQ*R 

AC«(X1-XR)/RC 

BC«(Yl-YR»/RC 

CC»(Zl-ZR)/RC 

C DETERMINES  WHETHER  RAY  STRIKES  ENCLOSED  SURFACE 
XSSaXP 
YSSaYP 
ZSS=ZP 

CALL  TRSCYU,XSSiYSSfZSS) 

XRRaXR 

YRRsYR 

ZRRaZR 

CALL  TKSCY( I tXRKfYRRfZRR) 

P1*XRR-XSS 

P2»YRR-YSS 

P3»ZRR-ZSS 

xiapxvns,!) 

Yi»PYvns,n 

zispzvns.i) 

CALL  TRSCYn,Xl,Yl,ZU 

Al*Xl-XSS 

Bl«Yl-YSS 
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Cl«21-ZSS 

IN«NVIIS) 

OC  1C  II=UIN 

iFni.EQ.iN)ic«: 
x2»pxvns,ic ) 

Y2»PYV(IS,IC) 

Z2«PZV(IS«IC) 

CALL  TRSCY(I*X2tY2tZ2) 

A2«X2-XSS 

B2-Y2-YSS 

C2«Z2-ZSS 

Q»Pl*(Bl*C2-B2*Cl)-P2*(Al*C2-Cl*A2)*P3*(Al*E2-ei#/2» 

IF(Q.GT.C.)RETLRN 

C RAY  STRIKES  SLRFACE  ON  ENCLOSED  SICE  OF  SURFACE  ECGE 
Al«az 
B1«R2 
Cl«C2 

1C  CONTINUE 

C RAY  STRIKES  ENCLOSED  SLRFACE 
ISK>I 
RETURN 
END 

/FORT 

C 

C 

C CONVERTS  CQCRCINATES  IN  CLINDRICAL  COORDINATES  INTO  RECTANGULAR  SPACE 
SUBROUTINE  TRSCY( I ,XRtYRtZR) 

COPMON/CYL/NCY*NSC(IO) *NSP ( 10 , 10 ) ,XC( 10 ) ,YC (10 ) , ZC ( 10 ) , AE ( 10 )»BE 1 1 
OOtCEdOltRCdO) 

CN«SQRm.>C£<n<**2) 

AN«*AE(nACE(n/CN 
BN«-BE(n*CEd)/CN 
Ap«Bem*CN-BN*cE(n 
6P»- ( AE  ( n*CN-AN*CE  d n 
CP«AEd)*8N-AN*BE(l) 

R0»AE(I)*(XR-XC(1  n*BEdM(YR-YC(m*CE(I»A(ZR-ZC(m 

XX*XC(I)*AEm*RO 

YY»YC(I  MBEd  )<R0 

ZZ«ZC(n+CE(I)*RO 

R«SQRT((XR-XX)**2*(YR-YY>**2*(ZR-ZZ»**2> 

AA«(XR-XX)/R 

Bi«(YR-YY)/R 

CC»IZR-ZZ)/R 

A»AN*AAtRN*fi6+CN*CC 

IF(A.GT.I.)Asl. 

1F(A.LT.-1.  IA»-1. 

ANG»ACOS(A) 

Q*AP*AA*BP*eB+CP*CC 
IF(O.LT.C.)ANG<-AN'G 
XR>ANG*RCd  ) 

YR»RC 
ZR>R*RC( I ) 

RETURN 

ENC 

/FORT 

C 

C 

C CCHPUTES  INCIDENCE  ANGLEt  REFLECTANCE,  AND  TRANSPITTANCE 
C NATURAL  LIGHT,  ACOITICN  OF  POLARIZATION  COMPONENTS  IGNORED 
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SUBRCCTINE  CCPP(AKG,R1 ,TT) 

CORHON/LINE/ASteStCSiXStYS.ZSfACfBCtCC 
C XNt  INDEX  OF  REFRACTION,  TX,  INTERNAL  TRANSNITTANCE 
CATA  XN.TX/1.S..92/ 

A»-AS*AC-BS*BC-CS*CC 

IF(ABS(AI.GT.l.)A=l. 

ANC=ACOS(A) 

ANGP-  ASIN(S1N(ANG)/XN) 

CA»CCS(ANG) 

SA>SIN(ANG) 

S1=SQRT(XN*A2-SA*A2> 

RO»(  ((CA-<H/(CA+Sin*A2*nCA*(XN**2»-Sl»/ICA*(XN**2»*SinAA2»/2. 
TO*( l.-RCI*CA/CCS(ANGP» 

CA>CCS(ANGP) 

SA>SIN(ANGP) 

Sl»SORT(XN**2-SA**2l 

RI«( ((CA-S1I/(CA*S1) »**2^( (CAA(XN*A2)-Sl)/(CAA(XN**2)4Sin**2>/2. 
TI»(l.-RI)*CA/CCS(ANG) 

TT«TC*TI*TX/a.-<RIATXIA»2> 

RT»RC*RI*TT 

RETLRN 

ENC 

/FORT 

C 

C 

C CCHPLTES  ARCSIN 

FUNCTION  ASINIXI 
Y«X 

AX-ABS(Y) 

IF(AX«GE.1.CI  GO  TO  A 

AC«ATAN(Y/SGRT(l.-Y*Yn 

ASIN*AC 

RETLRN 

A IF(AX.GE.l.CCCl)  GO  TO  10 
ASIN»l. 57079 
IF(Y.LT,C.)ASIN»-ASIN 
RETURN 

IG  NRITEdtlUX 
STOP 

11  FCRHATIBHERRCR  * ,23HARCSIN  ARCLPENT  .GT.l.  ,6FANC  > ,E16.8) 

END 

/FORT 

C 

C 

FUNCTION  ACCSIXI 
Y»X 

AXsABSIYt 

IFIAX.GE.l.CIGCTO  A 
IF(Y.EQ.C.)GCTC  3 
AC»ATAN<  SQRT(1.C-Y*Y»/Y) 

1 IF(Y.LT.C.  )GCTC  2 
ACCS=AC 

RETURN 

2 ACCS=AC*3.1A15926 
RETURN 

3 ACCS»1.57C79t3 
RETURN 

A IF(AX.6E.1.CC0CI )GOTO  10 
AC«C. 

GOTO  1 


ARCCCS  2 
ARCCCS  3 
ARCCCS  A 
ARCCOS  5 
ARCCCS  t 
ARCCOS  7 

ARCCOS  9 

ARCCCSll 

ARCCCS13 

ARCCOSIA 

ARCC0S15 

ARCCCSK 
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10  hRITEdtinx 
STOP 

11  FCRPAT(8HERRCR  * .23HARCC0S  ARGIMEKT  .GT.l.  t6HANC  « fE16. 
ENC 

/FORT 


ORAhS  GRAPHIC  PICHRE  OF  CANOPY  IN  3 FOLD  LAYOUT 
SUBROUTINE  CRNCN 
COHMON/GAREA/1CFIL(6COO) 

COPHON/C AN/NT tNB f NCtNPtNA,KCtNV 1100 ItPXVI 100. 8 ),PYV( 100,8), 
Q,8) 

CONPCN/NORP/AXMICO)  ,AYN(1CC)  .AlNIlOO) 

COHMON/PILCT/XC,YC,ZO 
COPMON/FACT/SX 
SX>2.122 
PRITE(1,SS8) 

898  F0RPATI2X, 'CONSOLE  NO  1,  TURK  OK') 

PAUSE 

CALL  GINUCCC) 

IX«(Y0-S.!)*SX41OC« 

IV»(20-111.52)ASX+1C0. 

CALL  GBEG(1,IX,IY) 

CALL  EPARK 

IY«(X0^10A. )*SX«433.AA 
CALL  GCPV(2,1,IX,IV) 

1X>I2G-111.S2)*SX«680.A2 
CALL  GBEG(3,IX,IY) 

CALL  ETtC 
NE*3 

CO  10  I>1,N0 

C ENTITY  IS  CANOPY  SURFACE— FENCE  OR  TRANSPARANT 
QS«-AXN(I ) 

QF»-AYN(I ) 

QT>4AZN<1) 

IX«(PYV(1,1)-5.5)*SX*1C0. 

IFm«LE.NC*ANO.QS.LT.O*).OR.n«6T.KC.ANG.CS*GT.O.))GOTO  2 
C SURFACE  FACES  VIENER  FROP  SIDE  VIEP 
IY*(PZV(I,I)-1U.S2)*SX«1C0« 

NE>NE«1 

CALL  GBEG(NE,IX,IY) 

IFII.LE. NOCALL  GPUT ( 3 ,130 ,1 ,2) 

IF(I.EQ.NA)CALL  GPUT(3,130,l,2) 

CALL  LINSU) 

2 CONTINUE 

IY«(PXVn,l)4lOA.  )*SX1A33.4A 

IF((I.LE.NC.ANO.QT.LT.O.).OR.(I.GT.KC.AND.QT.GT.O.))GOTO  3 
C SURFACE  FACES  VIENER  FROP  TOP  VIEN 
NE«NE«1 

CALL  GBEG(NE,IX,IY) 

IFII.LE. NOCALL  GPUTt3,l30,l,2) 

IF! I.EQ.NA)CALL  GPUT(3,130,l,2) 

CALL  LINTII) 

3 CONTINUE 

IFII  I>.LE.NC.ANO.QF.LT.O.).OR.<I.GT.KC.ANO.CF.GT.O.) )GOTO  10 
C SURFACE  FACES  VIENER  FROP  FRONT  VIEN 
IX><PZVII,l)-111.52)*SXt690.42 
NE*NE«1 

CALL  6BEG(NE,IX,IY) 

IFII.LE. NOCALL  GPUT(3,130,1,2) 


ARCC0S18 
) ARC00S19 


PZVIlOO 
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IF(I.EQ.NA)CAU  GPUT(3»13C»1 ,2) 

CALL  LINFd) 

1C  CONTINUE 
CO  20  I=>1»NE 
CALL  GFCNM) 

20  CONTINUE 
CALL  GSTART 
RETURN 
ENC 

/FORT 

C 

C 

C NARK  EYE  POSITION  IN  CANOPY 
SUBROUTINE  EPARK 
CALL  GPUT(6t50,-lC,-5» 

CALL  GPUT(7t70»0»10) 

CALL  GPUT(a,50»10«-5) 

CALL  GPUT(9,70f-7,-l0) 

CALL  GPUT(lCt50«0»20) 

RETURN 

ENO 

/FORT 

C 

C 

SUBROUTINE  ETIC 
CALL  GPUT(6t70f0«*S) 

CALL  GPUT(7tS0»0tI0l 
CALL  GPUT(8«70«-5t>5) 

CALL  GPUTdtSOflOiO) 

RETURN 

ENC 

/FORT 

C 

C 

SUBROUTINE  LINSU) 

CONNON/CAN/NTtNBtNCfNP,NA*NCtNV(IOO)»PXVaOO(8)tPYV(lOO»8)tPZVnOO 

Qt8) 

CONNON/FACT/SX 

CALL  GPUTOf I3Cf2tOI 

NK>NVm 

GO  10  K«1,NK 

K1>K«1 

IF(K.E.^.NK)  Kl«l 
IX«IPY'MItKn*PYV(lfK))«SX 
lY><PZV(I«KU-P2V(ItKn«SI 
CALL  GPUT(K«5t53»IX«IV) 

10  CONTINUE 
RETURN 
ENC 

/FORT 

C 

C 

SU8R0UT1NE  LINFd  I 

C0NM0N/CAN/NTtNB«NCtNP»NA,NDfNV(IOO)tPXVdOOi8)fPYV(lOO»8)«PZV(lOO 

Qf8l 

CONNON/FACT/SX 

CALL  GPUTI3fl30f2,0l 

NK«NVd) 

CO  10  K>1,NK 
K1>K*1 


5 

i 
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IF(K.EQ.NK)  Kl«l 
lX>(PZV(I»Kl)-P2V(ItKH*S)l 
IY-<PXVfI,Kl)-PXV(I tK))*SX 
CALL  GPUT(K«St53,IXt;Y) 

10  CONTINUE 
RETURN 
END 

/FORT 

C 

C 

SUBROUTINE  LINT(I) 

COMHON/CAN/NT*NB*NC«NP*NA«KO«NVaOO)»PXVaOOt8)  tPYVa00*8)tPZVfl00 

Qt6) 

CONPON/FACT/SX 

CALL  GPUT(3tl30«2«OI 

NK«NVm 

00  10  K>1,NK 

K1«K*I 

IFIK.EQ.NKI  Kl«l 
ix«(PYV(itxn*-PYV(itxn*sx 
lY«(PXVn«Kl)>PXV(I*KH*SX 
CALL  GPUT(K4S,S3«IXtIY) 

10  CONTINUE 
RETURN 
ENC 

/FORT 

C 

C 

C ORAUS  PERFECTIVE  OF  CANOPY  AKC  ENTRY  ANC  REFLECTION  POINTS  PILOT  EYE 
C POSITION 

SUBROUTINE  PERP(IES) 

COMNON/GAREA/IOFIL(6000) 

CONN0N/PER/N0tNS(S0l(PXS(SO»20)»PYStSOt2O)»PZS(50»20) 

DATA  AX«BXfCX/10,*S.t912./ 

CALL  TRSOC 
CALL  GINUCCO) 

1XE«CX 

IYE«CX 

CALL  GBEG(ltlXEflYE) 

CALL  ETIC 
CALL  GEONU) 

NE»2 

IK>10 

CO  10  I*1«N0 
NK*NSUI 

IFIlK.GT.e)  CALL  GBEG  INE *IXE »I YE ) 

1K«6 

00  9 K«1.NK 
Kl«K4l 

IFIK.EQ.NK)  Kl«l 
C CHECK  EDGES  FCR  HIDDEN  LINES 
CALL  CLIPLdK.ltK.Kl) 

5 CONTINUE 

IFIIK.EQ.6)  GO  TO  10 
CALL  GEON(NE) 

NE-NE41 
10  CONTINUE 

IF(IES.EQ.C)  GO  TO  21 
CALL  G8EG(NE4l,IXE>lYE) 

CALL  GPUT(5,176C»Ci. FALSE.) 
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IFIIES.E0.2)  CALL  GPtT(Stl750*0»0l 
CALL  POINT(IES) 

CALL  GECNIKE^l) 

21  CALL  GSTART 
RETURN 
END 

/FORT 


CONVERTS  OBJECT  COORDINATES  TO  DISPLAY  COORDINATES  » REMOVES  LINES  BEHIND 
VIEkERf  SURFACES  QLTSIOE -VIEWING  BOX,  SURFACES  NOT  FACING  VIEWER 
SUBROUTINE  TRSDC 

CONMON/CAN/NTtNB«NCtKPtNA«NCO*NV(100)tPXVI100t8)fPYV(100t8)fPZV(10 
QC«8) 

CONNON/NORN/AXNUCO)  tAYNUCO)  *AZNC100) 

CONNON/PILCT/XOtYCtZO 
C0NH0N/LINE/Xl«YltZltX2«Y2*22*X3fY3,Z3 
COMMON/PER/NO, NS(SO) tPXS(50 t20) *PYS I$0t20) t PZS(S0t20) 

DIMENSION  XS(lC0t8)yYS(100t8)(ZS(100»8) 

DATA  A«B/10.t5./ 

N0>1 

NSS*NB-»1 
CO  10  I*N!StND0 
IF(I«EQ.NA)GCTC  1C 

Q«AXNm«(pxv(i«ii-xo)^AYNm«(PYvn»i)-YO)4AZNm»(pzvntn-zo) 
IFIO.GT.O.IGCTO  10 
C SURFACE  FACES  VIEWER 
KN«NV<I) 

DO  2 K«1,WN 

CALL  PLAC(PXVn«K)tPYV(ItR).>PZVn»K|,XS(NOyK),YS(NO(K)>ZSINOyRII 
C CONVERTED  TO  VIEWER  COORDINAHS 

2 CONTINUE 
K0«0 
IXP«0 
1XN«0 
IZP«0 
IZN*0 

DO  5 K*1»KN 
Kl»K4l 

IFIK.EQ.KNI  Kl«l 

IF(YS(NO«K).LT,O..ANO.YS(RO»K1).LT.O.)  GO  TO  5 
IF(YS(NO«K).LT.O..OR.YS(NOfKl).LT.O.)  GO  TO  3 
XlaXSINOtK) 

Y1>YS(N0»KI 

Zi«ZS(NO,K) 

Q>1 

GO  TO  4 

3 CONTINUE 

CALL  INTEPIXS(NO«K)rYS(NO«K)yZS(ND»K)»XS(NC*Kl)tYS(NCyKlUZSINO»Kl 
ClyQI 

C SURFACE  EDGE  EXTENDING  BEHIND  VIEWER  TRIMMED 

4 KD*KD4l 

PXS(ND»K0)«X1*A/(Y1*BI 

PYS(ND,KDI»-1./Y1 

PZS(ND»K0)«Z1*A/(V1«B) 

CALL  TEST (PXSINDt KOI yPZS (NO .KOI ,IXP,IXN»IZP* IZN I 
C POSITIVE  POINT  OR  CROSSING  LINE  FROM  BEHIND 
IFiQ.GT.O.)  GO  TO  5 
K0«K041 

PXS(ND»KDI>X2*A/(Y2*B) 
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PYS(N0,KD)«-1./Y2 

PZS(N0tK0)«22*A/<Y2*ei 

C LEAC  POINT  FOR  LINE  CROSSING  PROP  IN  FRONT 

CALL  TEST(P)(S(ND«KO)  tPZS(KDtKO)tIXP,IXN»IZP»IZN) 

C VERTEX  CONVERTED  TO  DISPLAY  VIEkBOX  COORDINATES 
5 CONTINUE 

IF(KO.EQ.C)  GO  TO  10 
C SURFACE  NOT  BEHIND  VIEWER 

IFIIXP.GE.KO.OR.IXN.GE.KO.OR.IZP.GE.KO.OR.IZN.GE.KC}  GO  TO  10 
C SURFACE  PARTIALLY  OR  COHPLETELY  WITHIN  VIEW  BOX 
NS(ND)*KO 
ND«N041 
10  CONTINUE 
N0*N0-1 
RETURN 
END 

/FORT 

C 

C 

SUBROUTINE  PCINT(IES) 

DATA  AX«BX«CX/10,«S.,$12./ 

CALL  GET0EVILUN«*FINKAHS10) 

CALL  SEEK(LLN«C) 

REAOlLUNINVkC 
CALL  SEEK(LUN«NVNC> 

REACILUNINVW 

IK>6 

DO  20  I«ltNVW 

REA0<LUN)AlV»A2VfASPfBSP*CSP*IPtXPtYP«ZP,AIPtBIPtCIPtARP,CRP,CRP« 
QR«IS«XRfYR,ZR«AItBI*CItARtBR»CRtT 
IFIIES.EQ.l)  GO  TO  30 
60  TO  40 
20  CONTINUE 

CALL  OEVTILlNtOyOi 
RETURN 

C ENTRENCE  POINTS  FOR  EXTERNAL  RAYS 
30  CONTINUE 

CALL  PLACIXRyYRyZRfPXyPYyPZ) 

IFIPY.LE.C.)  CO  TO  20 

XX>PXPAX/(PY«8Xt 

Z2»PZ*AX/IPY«BX» 

IFIABS<XX).GT,l..CR.ABS(Z7).GT.l.)  GO  TO  20 

IXaXX«CX*CX 

JY.aZZACXKX 

GALL  6PUTnK,43»IX,IY) 

IK«IK«1 
GO  TO  20 

C PRIMARY  REFLECTION  POINTS 
40  CONTINUE 

CALL  PLACIXPtYPtZPyPXyPYyPZ) 

IFIPY.LE.C.)  GO  TO  20 
XX«PXAAX/(PY*BX) 

ZZ*PZ«AX/(PY<»BX) 

IF|ABSIXX).GT.1..0R.ABStZZ).GT.l.)  GO  TO  20 

IX*XX*CX4CX 

IV»ZZ*CX*CX 

CALL  GPUT(IK(1G0»IX>4«0) 

CALL  GPUT(lK*ltllC,lY-4,0) 

IC»-ALOGlC(T) 

IC»IC4l76 
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CALL  GP(JTnK-»2t90tIC«0> 

1K>IK*3 
6C  TO  20 
ENC 

/FORT 

C 

C 

C DETERMINES  VISIBLE  PORTION  OF  EDGE  BY  REMOVING  HICCEN  LINE 
C ALL  CONSTRAINT  SURFACES  CONVEX 

SUBROUTINE  CLIPL(IK«IV,KV«KV1) 
COMMON/LINE/AV«BVtCV*XVtYVtZV»ACfBC»CC 
C0MM0N/P£R/NC«NS(50) ,PXS (50 *20 ) »PYS (50,20) »PZS ( 50 »20 ) 
COMMON/DRA(</NSC*RV*ROVaO)*RFV(10) 

DATA  CX/512./ 

XV«PXS(IV«KV) 

YVPYSnVtXV) 

ZV*PZS(IV«KV) 

RV«SQRT(  (PXS(IV»KVII-XV»**2*(PYSIIV,KV11-YVMA24(FZS(IV,KV1I-ZV|M* 
02) 

AV«(PXS(IV*KV1)-XV)/RV 
BV*(PYS(IV«KV1)-YV)/RV 
CV*(PZS< IV«KV1)-ZV)/RV 
R0>0. 

RFsRV 

CALL  CHKS(-1.,0.«C.,R0»RF) 

CALL  CHKS(l.tC.fO.«RO«RF’/ 

CALL  CHKS(C.,0.tl.»RO*RF) 

CALL  CHKS(0.t0.«-l.«R0»RF) 

IF(RF.LT.RC)  RETURN 
C ECGE  CNECKEO  AGAINST  VIEhBOX 
NSC«1 
R0V(1)«R0 
RFVm»RF 
CO  10  I«ifNO 
IF(tv.EQ.I)  GO  TO  10 

C FIND  APPARANT  INTERSECTION  POINT  OF  EDGES 
ISC*0 
KN«NSm 
00  5 K>lfKN 
Kl«K+l 

IF(K.EQ.KN)  Kl«l 
XE>PXS(ItK) 

YE«PYS(IfX) 

ZE«PZS(I«K) 

RE«SQRT( (PXS(ItKl)*XE)«*24(PYS(IfKl)>YE)*«2*(PZS(I»Kl)>ZE)**2) 
AE«(PXS(ItKl)-XE)/RE 
BE-IPYS(IfKl)-YE)/RE 
CE«(PZS(IfKl)-ZE)/RE 
Q»AE*CV-AV*CE 
IF(Q.EQ.C.)  GO  TO  5 
C EDGES  NOT  PARALLEL 

RIV»(CE*(XV-XE)-AE*(ZV-ZE))/C 
IFIAE.EQ.C.)  GO  TO  2 
RIE*IXV-XE«AV*RIV)/AE 
GO  TO  3 

2 CONTINUE 

RIE»(ZV-ZE4CV*RIVI/CE 

3 CONTINUE 
YIE=YE*BE*RIE 
YIV»YV*BV*RIV 


50 


uuuuu 


IF(YI£.6T.Y1V)GCTC  5 
C CONSTRAINT  cOGE  IN  FRONT  OF  TEST  ECGE 

IF(RIE.GT«O..AND.RI£.LT.RE)  GO  TO  7 
5 CONTINUE 

IFIISC.lt. GO  TO  10 
IFIRO.lt. 0.. AND. RF.GT.RV)  RETLRN 
C TEST  EDGE  NOT  FICOEN  BEHIND  CONSTRAINT  SURFACE 
IFIRO.GT.RV.OR.RF.LT.O.)  GO  TO  10 
C TEST  EDGE  PARTIALLY  BLOCKED  BY  CONSTRAINT  SURFACE 
CALL  CKLINIRCtRF) 

IFINSC.EQ.O)  RETURN 
GO  TO  10 
7 CONTINUE 

C CONSTRAINT  EDGE  LOCATEO 
ISC-1SC«1 

IF(lSC.EQ.l)  R1*RIY 

IFIISC.EQ.2I  CALL  OROLNIRl »RIVtRO»RF) 

CO  TO  S 
10  CONTINUE 

C VISIBLE  PORTION  OF  EDGE  REPAINS 
CO  20  I>1«NSC 
R0>R0V(1I 
RF-RFVm 
XS«XV4R0*AY 
z$>z;v«ROACv 
XF«XV4RF*AV 
2F«2V4RF*CV 
IXS«IXS4l.l*CX 
IZS«I2S4l.)*CX 
1XF«IXF4I.IACX 
IZF«(ZF«1.)ACX 
CALL  GPUTUNtlCCflXS.O) 

CALL  GPUTIIK4iaiOfIZSfOI 

IDXalXF'IXS 

10Z«1ZF*IZS 

CALI.  GPUTnK42tS3tI0X«I0Z) 

1K»IK43 
20  CONTINUE 
RETURN 
ENC 

/FORT 


CONVERTS  OBJECT  SPACE  COORDINATES  INTO  EYE  SPACE  COORDINATES 
PILOT  EYE  OIRECTEB  IN  Y-Z  PLANE  OF  AIRCRAFT  COORCINATES  AND  5-DECREES 
ABOVE  NEGATIVE  Y-AXIS 

SUBROUTINE  PLAC (PXV»PW»PZV tPX tPYtPZ I 
COPNONA IlCT/XOfYO(ZO 
CONMON/ANG/AKfBN 
DATA  PI/3.14159/ 

ANl«AN*Pi/180. 

BNl-BNPPI/IEO. 

PX»(PXV-XCI*SINCBNll-iPYV-YC»*COS(BNll 

PY-  (PXV-XCI«CCS(BNll<»SIMANl)-iPYV-VOI*SIN(eNU*SlNIANll-(PZV-ZOI 
Q4C0SIAN1I 

PZ»-IPXV-YOI*COS«BN1)*COSIAM)-CPYV-YOIPS1NIBM»ACOS(AN1MIPZV-ZOI 
Q*SIN(AN1I 
RETURN 
ENC 

/FORT 
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c 

C TEST  SUFACE  VERTEX  FOR  POSITIOK  OUTSICE  VIEN80X 
SUBROUTINE  TEST (XyZ»lXPtI'XKtXZP*IZM 
IFIX.GT.l.)  IXP«IXP«1 
EF(X.LT.-1.)IXN>1XN«1 
IFIZ.GT.l.)  IZP«IZP*1 
IF(Z.LT.-1.)IZN>1ZN>1 
RETURN 
ENC 

/FORT 

C 

C 

C CALCULATES  INTERSECTION  POINT  FOR  LINE  WITH  PLANE  NORPAL  TO  Y-AXIS 
SUBROUTINE  INTEPf Xl*Yl,Zl^X2tY2 »Z2»B) 
COPNON/Ll^E/XStYS«ZStXF»YFtZF»XC*YCfZC 
CATA  YO/1./ 

R«SQRTM  X1-X2  ( YI-Y2 1♦♦2♦  <Z1-Z2  )**2) 

A«(X2-X1)/R 
B-IY2-YII/R 
C«(Z2-Z1I/R 
IFIB.EQ.O.)  RETURN 
C ECGE  NOT  PARALLEL  TO  PLANE 
YI«YO 

RI*(YI-Yl)/8 
XI>X1«A*RI 
ZI«Zl*C*RI 
IFIB.LT.C.I  GO  TO  3 
C LKNE  ORIGINATES  BEHINC  VIEMER 
XS>XI 
YS-YI 
ZS>ZI 
XF>X2 
YF»Y2 
ZF«Z2 
RETURN 
3 CONTINUE 

C LINE  ORIGINATES  IN  FRONT  OF  VIEWER 
XS«X1 
YS«Yl 
ZS>Z1 
XF>XI 
YF«YI 
ZF-ZI 
RETURN 
ENC 

/FORT 

C 

C 

C CHECKS  EDGE  AGAINST  VIEWING  BOX  SIDE 
SUBROUTINE  CHKS (AN»BN*CNfROfRF) 
COMMON/LINE/AVtBVtCVfXVyYVfZVyACfBCyCC 
XN«-AN 
YN«-BN 
ZN>-CN 

Q«AV«AN«BV*8N-»CV*CN 
IFIO.EQ.C.)  GO  TO  5 

RI»HXN-XV)*AN*(YN-YV»*BN«IZN-ZV»*CM/Q 
IFIQ.LE.O.)  GO  TO  3 

C ECGE  DIRECTED  INTO  VIEWBOX  FRCP  OUTSIDE 


52 


O O O ^ V)  u u 


IF(RI.GT.RC)  RC«RI 
RETURN 
3 CONTINUE 

C EDGE  DIRECTED  OUT  OF  VIEMBOX  FROF  INSIDE 
IFIRI.lt. RF)  RF«RI 
RETURN 
5 CONTINUE 

C EDGE  PARALLEL  TC  SIDE 

IF(ABS(XV).GT.1..0R.ABS(ZV).GT.1.I  RF«RO-l. 

RETURN 

END 

/FORT 


ORDER  EDGE  INTERSECTION  POINTS  ACCORDING  TO  LOM  AND  FIGH  VALUES 
SUBROUTINE  ORDLNI R1 tR2 »RO«RF ) 

RO-Rl 


RF«R2 

IFIR1.lt. R2I  RETURN 

R0«R2 

RF>Rl 

RETURN 

ENC 


FORT 


CHECKED  PARTIALLY  BLOCKED  LINE  FOR  VISIBLE  SE6PENTS 
SUBROUTINE  CKLIN(R1,R2) 
COMNON/ORAM/NSCtRVtROVaO)tRFV(10) 

DIMENSION  ROaO)tRF(lO) 

IF(R1.GT.0..AND.R2.LT.RV)  GO  TO  20 
C VERTEX  OF  TEST  EDGE  BEHIND  CONSTRAINT  SURFACE 
IF(R2.6T«RV)  60  TO  5 
C LOM  END  OF  EDGE  HIDDEN 
RV0-R2 
RVF>RV 
GO  TO  10 

C HIGH  END  OF  EDGE  HIDDEN 
5 CONTINUE 
RVC«0. 

RVF«R1 
10  CONTINUE 

DO  12  I«ltNSC 

Rom«Rovm 

RFm-RFvm 

IF(RFVm.LE.RVO)  RFID*-!. 

IFIROVm.LT.RVOi  ROm>RVO 
IFIRFVm.GT.RVFI  RF(I)>R>VF 
IFIROVm.GT.RVF)  RFm>-l. 

12  CONTINUE 
GO  TO  30 

C CONSTRAINT  SURFACE  SEPARATES  EDGE  INTO  TNO  VISIBLE  ENDS 
20  CONTINUE 
RVC«R1 
RVF«R2 
K«0 

CO  22  I«ltNSC 
K«K«1 

RO(K)«ROVm 

RF(K)«RFVm 
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IFIRFVm.GT.RVO.ANO.RFV(I).LT.RVF)  RF(Kl«RVO 
IFiROVn  ».GT.RVO.ANO.RFV(l).LT.RVF)  RF|K)«-1. 
IF(RFVn).G1.RVF.ANO.ROV(I).LT.RVO)  60  TO  27 

IF(RFVn).GT.RVF.AND.R0Vm.6T.RV0.AK0.R0V(n.LT.RVF)  RO(K)«RVF 
22  CONTINUE 
NSC«K 
GO  TO  30 
27  CONTINUE 
RFIK)«RVC 
K«K*l 
RO(K)*RVF 
RF(K)>RFVm 
GO  TO  22 

C ARRANGE  REMAINING  VISIBLE  SEGMENTS 
30  CONTINUE 
NK«0 

00  32  I^lfNSC 
IFIRFIII.LE.C.)  GO  TO  32 
NK«NK4l 
ROV(NK)>RCm 
RFV(NK)>RF(n 
32  CONTINUE 
NSC«NK 
RETURN 
ENC 

/FORI 


SUBROUTINE  PERCdA.IB) 

C0NN0N/AN6/AN,BN 

AN«IA 

BN«IB 

CALL  REAOF 
CALL  PERP(l) 

PAUSE 

CALL  GS3UU 
PAUSE 

CALL  PERP(2) 

PAUSE 

CALL  GS3l(l) 

PAUSE 

RETURN 

END 


DRAWS  GRAPHIC  PICTURE  OF  CANOPY  LAYOUT  AND  ENTRYt  PRIMARY  REFLECTION  POINTS 
SUBROUTINE  ORWCP(lEflR) 

COMMON/FACT/SX 

CALL  0RNCFINEflXl»IYl(IY2»lX3) 

CALL  GBEG(NE«ltIXl,lYU 
CALL  GPUII5»176CfO». FALSE.) 

CALL  6BEG(NE>2tIXl.lYl) 

CALL  GPUT(5tl750»C«C) 

CALL  GREG(NE^3«IXI,1Y2) 

CALL  6PUTI5, 1760. Of. FALSE.) 

CALL  6BEG(NE«AfIXl.IY2) 

CALL  GPUTIS. 1750.0.0) 

CALL  6BE6(NE«5.IX3.IY2) 

CALL  GPUTIS. 1760*0.. FALSE.' 
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CALL  6BEG(NE«6*IX3»IY2) 

CALL  GPUT(5tl750tC,0) 

IKP*6 

IKR«6 

CALL  GETOEVILUNf 'FINKAH',10) 

CALL  SEEKILLNtC) 

READlLUNINVIiC 
CALL  SEEK(LLN*NVWO) 

REAOILUNINVb 
00  20  I*ltNVh 

REA0{LUmAlV»A2VtASPtBSP«eSP,IPtXP,YPfZPtAIP»eiP,CIP«ARP,eRPfCRPt 
QR.IStXRtYRtZRtAltBItCI tARtBRtCRfT 
IFUS.EQ.IE.ANO.IP.EQ.IRIGQTO  30 
20  CONTINUE 

CALL  0EVT(LUK«0,Ot 
PAUSE 

CALL  GS3K1) 

PAUSE 

RETURN 

C ENTRENCE  POINTS  FOR  EXTERNAL  RAYS  AND  REFLECTION  POINTS  TO  PILOT'S  EYE 
30  CONTINUE 
IC«-AL0G1C(TI 
IC«IC«176 
CALL  6ENTINE41) 

IXR«fVR-3C.I*SX 

IVR>(ZR-1CA.A8I«SX 

CALL  GPUTdKRfAStlXRtlYR) 

IXP«<VP-3A.I«SX 
IVP«(ZP-1C6.AB)«SX 
CALL  6ENTINE42) 

CALL  GPUT<IKP*100.IXP-A«0) 

CALL  6PUTnKP«ltll0fIYP*4t0) 

CALL  6PUTUKP«2t9CtIC»0) 

CALL  GENTINE43) 

IVR«(XR4121.44I*SX 
CALL  6PUTnKR«43*IXRtIYR) 

IVP«(XP4121.44)4SX 
CALL  GENT(NE44) 

CALL  GPUTnKP,l00tIXP-4t0) 

CALL  6PUTIIKP4l,110fIYP-4«0) 

CALL  GPUTnKP42,90«IC*0) 

IXR«(ZR4IS.Sn*SX 
CALL  GENTINE4S) 

CALL  6PUTnKR«43tlXRtIYRI 

IXP«tZP4l3.9I)«SX 

CALL  GENTINE4A) 

CALL  6PUTnKP«100«IXP-4*0) 

CALL  6PUTnKP4l,110tIYP-4,0) 

CALL  6PUTnKP42t90fIC»0l. 

IKR>IKR41 
IKP«IKP43 
60  TO  20 
ENO 

/FORT 


ORAMS  GRAPHIC  PICTURE  OF  CANOPY  LAYOUT  AND  ENTRY  OR  FRIPARY  REFLECTION  POINTS 
ON  ALL  SURFACES 

SUBROUTINE  ORNCT(IE) 

CONMON/FACT/SX 


CALL  0RUCF(NE»lXl»m«IY2,I)l3) 

CALL  GBEGINE4l»IXl,IYl) 

CALL  GPUT(S»176CfC».FALSE.) 

I^(IE.E().2)CALL  GPUT(S*17SO»0»0) 

CALL  GBEG(NE«2>IXltIV2l 
CALL  GPU1(5tl76C,0». FALSE.) 

IF(IE.EQ.2)CALL  GPLT(5»175Ct0»0) 

CALL  GBEG(NE«3.IX3tIY2) 

CALL  GPU.T(5«17£0*C». FALSE.) 

IFnE.EQ.2)CALL  GPUT (5*1750*0 »0) 

IKP>6 

IKR>6 

CALL  GET0EV(Lt'N**FINKAH«*10) 

CALL  SEEK(LLN*0) 

REAC(LUN)KVtiC 
CALL  SEEK(LLN,NVMC) 

REA0(LUN)KVIi 
00  20  I>l,NVh 

REA0<LUN)A1V*A2V*ASP*BSP*CSP*IP*XP*YP*ZP*AIP*BIP*CIP*ARP*BRP*CRP* 
QR«1S*XR*YR*2R«AI*BI*CI*AR*BP*CR*T 
IF(IE.EQ.2)G0T0  3C 
60  TO  40 
20  CONTINUE 

CALL  0EVT(LUN*0*0) 

PAUSE 

CAU  GS3im 

P«rUSE 

RETURN 

C PRINARY  REFLECTION  POINTS  TO  PILOT’S  EYE 
30  CONTINUE 

1C«-ALOG10(T) 

IC«IC*176 
CALL  GENTINE^l) 

IXP»(YP-36.)*SX 

IYP«(ZP>IC6.48)*SX 

CALL  GPUT(IKP*100*IXP>4*0) 

CALL  GPUT(IKP4l*ll0*IYP-4*0) 

CALL  GPUTnKP«2*9C*IC*0) 

CALL  GENT(NE«2) 

IYP«(XP*121.44)«SX 

CALL  GPUT(1KP*100*IXP-4*0) 

CALL  GPUTIlKP+l,ll0*IYP-4*0) 

CALL  GPUT(IKP»2.9C*IC*0) 

IXP*(ZP4l5.Sl)<»SX 
CALL  GENT(NE*3) 

CALL  GPUT(IKP*100«IXP-4*0) 

CALL  GPUT(IKP^l*ll0*lYP-4*0) 

CALL  GPUT(IKP«2*9C*IC*0) 

IKP>IKP«3 
GO  TO  20 

C ENTRENCE  POINTS  FOR  EXTERNAL  RAYS 
40  CONTINUE 

IXR«(YR-36.)*SX 
IYR«(ZR-1C6.48)*SX 
CALL  GENT(NE>1) 

CALL  GPUT(IKR*43*1XR*IYR) 

IYR*(XR«121.44)*SX 
CALL  GENT(NE«2) 

CALL  GPUT(IKR*43*IXR*IYR) 

IXRa(ZR*l5.9l)*SX 
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o o o 


CALL  GENT(NE+3) 

CALL  GPUT(IKR,A3tIXR»lYR^ 

IKR=IKR«1 

GO  TO  20 

ENC 

/FORT 


DRAWS  GRAPHIC  PICTL'RE  OF  CANOPY  FPAPE  LAYOUT 
SUBROUTINE  ORIiCF(NEtIXl,IYl  »I  Y2  *1X3) 

COMNON/GAREA/IOFIL(8COO) 

COHMON/CAN/NT»NB*NC«NP*NA*ND»KV(100)«PXV(100»8)*PYV(100.8)*PZV(100 
0*8) 

C0HM0N/N0RR/AXN(100)*AYNI100)*AZN(100) 

COMMON/P ILOT/XO*YC«ZO 
COMMON/FACT/SX 
SXaA.86 
CALL  REAOF 
CALL  GINieCCC) 

1X1*(Y0>36.)«SX 
IY1*(Z0>1C6.48)«SX 
CALL  GBEG(1*IX1»IY1) 

CALL  EMARK 
IY2>(X04121.44)*SX 
CALL  GCPY(2«1*IX1«IY2) 

IX3>IZ0«13.91)*SX 
CALL  GBEG(3*IX3«IY2) 

CALL  ETIC 
NE«3 
NS«NC«1 
00  10  I«NS*NO 
IF(I.LE.NB)GCTO  10 
C ENTITY  IS  CANCPY  SURFACE 
QS«*AXNm 
QF*-AYN(II 
QT«4AZNm 

IX*(PYV(1*1)-36.)«SX 
IFIQS.GT.C.IGOTO  2 

C SURFACE  FACES  VIEWER  FROP  SIDE  VIEW 
IY«(PZV(I»1)-106.A8)*SX 
NE>NE41 

CALL  GBEGINEtlXflY) 

IFU.EQ.NAICALL  GPUT<3*130»1*2) 

CALL  LINSm 

2 CONTINUE 

IY<(PXV(I*1)4121.AA)*SX 
IFIQT.GT.O.IGOTO  3 

C SURFACE  FACES  VIEWER  FROP  TOP  VIEW 
NE>NE4l 

CALL  6eGGINE*lX,lY) 

IFII.EQ.NAICALL  GPUT(3*130*1*2) 

CALL  LINT(I) 

3 CONTINUE 
IF(QF.GT.C.)GOTO  10 

C SURFACE  FACES  VIEWER  FROP  FRONT  VIEW 
lX>(PZVn*l)«lS.91)*SX 
NE-NE41 

CALL  GBEGINEyIXtIY) 

IFII.EQ.NAICALL  GPUT  13*130 *1 *2 ) 

CALL  LINFin 


57 


1C  CONTINUE 
RETURN 
ENC 

C AAH  CANOPY  DATA 


c 

NUMBER  CF 

VERTICES, 

SURFACES- 

—CONSTRAINT 

, FLAT,  AND  CYLINCRICAL 

1E6 

16 

56  58 

59 

6A 

c 

NC. 

VERTICES  PER  SURFACE 

A 

A 

A A 

A 

A 

A 

A 

A 

A A 

5 5 

A 

A 

6 

A 

A 

A A 

A 

A 

A 

A 

A 

7 7 

7 7 

A 

A 

6 

A 

A 

A A 

A 

A 

A 

5 

A 

A 5 

7 7 

7 

7 

6 

A 

A 

A A 

A 

A 

A 

A 

A 

A 8 

6 6 

6 

6 

6 

c 

VERTICE  ASSIGNED  TC 

EACH 

SURFACE 

es 

90 

66  121 

125 

129 

133 

135 

139 

1A3  1A5 

IA9  153 

157 

137 

165 

1 

2 

3 A 

7 

1 

5 

12 

11 

21  15 

28  22 

32 

29 

33 

110 

35 

36  111 

38 

39 

38 

109 

AO 

A1  112 

53  A7 

60 

56 

61 

61 

89 

9A  S3 

ICO 

99 

98 

87 

77 

81  101 

68  71 

9A 

ICO 

85 

117 

llA 

115  122 

126 

130 

13A 

136 

lAO 

lAA  1A6 

150  ISA 

158 

136 

168 

2 

3 

A 5 

lA 

7 

A 

lA 

12 

20  16 

27  23 

31 

30 

36 

1C9 

36 

37  112 

39 

38 

111 

111 

A1 

110  no 

52  A8 

59 

55 

66 

62 

SO 

89  9A 

95 

ICO 

99 

98 

78 

82  102 

69  72 

93 

99 

86 

llA 

115 

119  123 

127 

131 

130 

137 

lAl 

lAO  1A7 

151  166 

159 

162 

163 

9 

10 

11  12 

12 

5 

3 

8 

9 

19  17 

26  2A 

30 

31 

109 

35 

A3 

AA  AA 

A6 

3A 

109 

37 

A5 

112  A2 

51  A9 

58 

56 

63 

63 

66 

65  70 

71 

76 

75 

7A 

79 

83  103 

70  73 

92 

98 

87 

9C 

66 

65  12A 

128 

132 

129 

138 

1A2 

139  1A8 

152  155 

160 

161 

166 

e 

9 

10  11 

5 

2 

2 

9 

10 

18  18 

25  25 

29 

32 

lie 

A2 

A2 

A3  37 

A5 

33 

3A 

36 

A6 

A5  A3 

50  50 

57 

57 

62 

6A 

65 

70  69 

76 

75 

7A 

88 

80 

8A  lOA 

65  7A 

91 

97 

88 

0 

C 

0 0 

0 

0 

0 

0 

0 

0 C 

165  156 

0 

C 

C 

c 

0 

0 0 

0 

0 

0 

0 

0 

17  19 

2A  26 

0 

0 

61 

c 

c 

C 0 

0 

0 

0 

35 

0 

0 AA 

A9  51 

56 

58 

C 

G 

c 

C 0 

0 

0 

0 

0 

0 

0 105 

66  75 

90 

96 

C 

C 

0 

0 0 

0 

0 

0 

0 

0 

0 0 

0 0 

0 

C 

C 

0 

0 

0 0 

0 

0 

0 

0 

0 

16  20 

23  27 

0 

0 

60 

0 

0 

C 0 

0 

0 

0 

0 

0 

0 0 

A8  52 

55 

59 

C 

Q 

0 

0 0 

0 

0 

0 

0 

0 

0 106 

67  76 

89 

95 

0 

C 

0 

0 0 

0 

0 

0 

0 

0 

0 0 

0 0 

0 

0 

0 

C 

0 

0 0 

0 

0 

0 

0 

0 

15  21 

22  28 

0 

0 

C 

0 

0 

0 0 

0 

0 

0 

0 

0 

0 0 

A7  53 

5A 

60 

8 

0 

0 

0 0 

0 

0 

0 

0 

0 

0 107 

0 0 

0 

0 

C 

0 

0 

C 0 

0 

0 

0 

0 

0 

0 0 

0 0 

0 

0 

C 

0 

0 

C 0 

0 

0 

0 

0 

0 

0 0 

0 0 

0 

C 

8 

c 

0 

0 0 

0 

0 

0 

0 

0 

0 0 

0 0 

0 

C 

C 

0 

0 

0 0 

0 

0 

0 

0 

0 

0 108 

0 0 

0 

0 

C 

c 

VERTICE  X-POSITICN 

2 

.cc 

2.00 

2.00 

2.00 

2 

.00 

2.00 

2.00 

-2.00 

-2 

.oc 

-2.00 

-2.00 

-2.00 

-2 

.00 

-2.00 

19.06 

18.70 

lA 

.25 

12.90 

11.20 

1A.76 

17 

.A2 

-18.00 

-17.60 

-13.86 

-13 

.30 

-13. 2A 

15.96 

-17.96 

8 

.50 

-11.50 

-12, OA 

9.08 

19 

.63 

19.63 

11. A6 

11. A6 

11 

.A6 

19.63 

19.63 

-19.63 

-19 

.63 

-11. A6 

11. A6 

- 

11. A6 

-19 

.63 

-IS. 63 

19.16 

18.86 

lA 

.22 

13.10 

11. 2A 

1A.78 

17 

.50 

-17.90 

-17.60 

-16.00 

-13 

.30 

-13.28 

15.96 

-17.96 

8 

.56 

-11. A8 

-12. lA 

9.16 

23 

.C8 

20.92 

18.08 

13.72 

22 

.72 

23.28 

23.20 

13.60 

13 

.CC 

15.16 

16.36 

22. 8A 

16 

.72 

-16.72 

-1A.28 

16.28 

13 

.60 

-13.60 

-9.20 

9.20 

10 

.56 

-10.56 

-11.56 

11.56 

-23 

.C8 

-20.92 

18.08 

- 

13.72 

-22 

.72 

-23.28 

-23.20 

-13.60 

-13 

.00 

-15.16 

16.36 

- 

22. 8A 

19 

.63 

13.26 

11. A6 

-11.66 

-13 

.26 

-19.63 

12.00 

12.00 

13 

.26 

-11. A6 

13.26 

-11.66 

-23 

.75 

-A. 

A. 

23.75 

-A 

• 

-23.75 

A. 

23.75 
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24. 

36. 

36. 

24. 

-36. 

-24. 

-24. 

-36. 

26. 

104. 

104. 

26. 

26. 

104. 

67.6 

61. 

61. 

67.6 

104. 

-26. 

-26. 

-104. 

-104. 

-26. 

-61. 

-67.6 

-67.6 

-61. 

-46. 

-24. 

-22. 

-56. 

24. 

46. 

56. 

22. 

-10. 

10. 

10. 

-10. 

61. 

61. 

-61. 

-61. 

-56. 

56. 

C VERTICE  Y- 

POSITION 

57. 5C 

57.50 

68.37 

71.57 

58.98 

58.98 

60.30 

57.50 

57. 5C 

68.37 

71.57 

58.98 

58.98 

60.30 

84.66 

82.92 

83.46 

86.98 

98.18 

98.66 

96.24 

84.66 

82.92 

83.46 

86.58 

98.18 

98.66 

96.24 

87.48 

87.48 

97.56 

97.56 

98.45 

103.49 

100.92 

110.92 

115.61 

121.40 

121.50 

98.45 

103.49 

100.92 

110.92 

115.61 

121.40 

121.50 

144.26 

142.52 

142.54 

146.00 

157.52 

156.40 

156.32 

144.26 

142.52 

142.54 

146. OC 

157.52 

158.40 

156.32 

148.00 

148.00 

156.66 

156.66 

59.20 

59.20 

70.14 

108.42 

110.90 

69.78 

116.49 

113.25 

139. 6S 

156.97 

156.97 

140.09 

58.30 

58.30 

67.18 

67.19 

69.57 

69.57 

108.01 

108.01 

112.13 

112.13 

145.44 

145.44 

59.20 

59.20 

70.14 

108.42 

110.90 

69.78 

116.49 

113.25 

139.69 

156.97 

156.97 

140.09 

115.52 

115.52 

114.82 

114.82 

115.52 

115.52 

113.01 

113.01 

103.49 

103.49 

121.40 

121.40 

57.5 

5.5 

5.5 

57.5 

5.5 

57.5 

5.5 

57.5 

59.2 

65.2 

132.61 

158.6 

65.2 

59.2 

158.6 

132.61 

186.61 

191.61 

226.61 

236.61 

186.61 

191.61 

189.61 

189.61 

180.61 

180.61 

191.61 

186.61 

236.61 

226.61 

191.61 

186.61 

189.61 

189.61 

180.61 

180.61 

214.61 

214.61 

214.61 

214.61 

214.61 

214.61 

214.61 

214.61 

198.61 

198.61 

198.61 

198.61 

198.61 

224.61 

198.6 

224.61 

214.61 

214.61 

C VERTICE  2-POSniON 

131.81 

138.66 

146.27 

141.77 

133.97 

133.97 

132.83 

131.81 

138.66 

146.27 

141.77 

133.97 

133.97 

132.83 

128.40 

130.34 

145.26 

148.54 

148.08 

136.54 

127.96 

128.40 

130.34 

145.26 

148.54 

148.08 

136.54 

127.96 

118.61 

118.61 

156.00 

156.00 

129.20 

148.00 

157.58 

160.51 

160.88 

139.20 

129.20 

129.20 

148.00 

157.58 

160.51 

160.88 

139.20 

129.20 

147.56 

149.40 

164.26 

167.64 

167.68 

156.20 

147.54 

147.56 

149.40 

164.26 

167.64 

167.68 

156.20 

147.54 

137.86 

137.86 

175.68 

175.61 

129.05 

141.41 

154.80 

172.40 

150.68 

130.85 

146.36 

174.59 

179.28 

179.24 

176.40 

150.80 

143.52 

143.52 

154.40 

154.40 

156.72 

156.72 

175.15 

175.15 

177.93 

177.93 

182.12 

182.12 

129.05 

141.41 

154.80 

172.40 

150.68 

130.85 

146.36 

174.59 

179.28 

179.24 

176.00 

150.80 

152.54 

152.54 

159.03 

159.03 

152.54 

152.54 

175.79 

175.79 

148.00 

148.00 

139.20 

139.20 

144.2 

129.1 

129.1 

144.2 

122.1 

127.05 

122.1 

127.05 

124.12 

124.12 

124.12 

124.12 

124.12 

124.12 

124.12 

124.12 

142.52 

140.52 

140.52 

142.52 

132.12 

132.92 

132.52 

132.52 

111.52 

111.52 

140.52 

142.52 

142.52 

140.52 

132.92 

132.12 

132.52 

132.52 

111.52 

111.52 

147.92 

147.92 

178.52 

178.52 

147.92 

147.92 

178.52 

178.52 

181.52 

181.52 

221.52 

221.52 

111.52 

132.52 

111.52 

132.52 

157.92 

157.92 

C CYLINDRICAL  DATA 
3 

2 2 
60  61 

1 

62  63 

64 

-84.74 

135.89 

130.93 

.0 

.9805 

.1968 

108.98 

84.74 

135.89 

130.93 

.0 

.9805 

.1968 

108.98 

C. 

146.86 

55.37 

1. 

.0 

.0 

127.63 

C PILOT  EYE 

POSITION 
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-1.5  142.33  171.2 

C CC-PILCT  EYE  POSITION 

-1.5  83.04  152.14 


